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PREFACE 


In  January  1985  the  Director  of  the  Science  and  Technology  Directorate  of  the 
Strategic  Defense  Initiative  Organization  (SDIO)  asked  members  of  the  IDA  research  staff 
to  investigate  the  feasibility  of  developing  a  y-ray  laser.  The  staff  determined  what  work 
had  been  done,  who  was  currently  working  in  the  field,  and  what  work  should  be 
encouraged  or  supported.  A  workshop  was  convened  for  researchers  directly  involved 
both  in  gamma-ray  laser  work  and  in  ancillary  fields  such  as  nuclear  structure,  radiation 
propagation  in  crystals,  Mdssbauer  Effect,  and  opdcal  lasers.  Next,  an  in-house  study  was 
undertaken  to  clarify  critical  issues  concenting  the  various  pumping  schemes  proposed  at 
the  workshop  as  well  as  systems  questions  about  the  f-ny  laser  as  a  wmking  device. 

The  proceedings  of  the  workshop  were  published  in  the  form  of  a  report  to  the 
Innovative  Science  and  Technology  Office  (1ST)  of  the  SDIO.  The  work  completed  in 
1985  was  presented  in  IDA  Paper  P-2021. 

In  1986,  the  in-house  work  focus^  on  extending  the  data  base,  the  nature  of 
superradiance  in  the  y-ray  laser  context,  and  a  detailed  investigation  of  the  upconversion 
pumping  scheme.  A  discussion  of  nuclear  systematics,  investigations  of  electron-nuclear- 
driven  pumping,  and  lifetime  measurements  rounded  out  that  study.  The  results  of  the 
FY  1986  effort  were  also  presented  in  an  IDA  paper,  P-2004. 

In  1987,  the  IDA  research  staff  looked  at  the  state  of  the  an  and  assessed  the 
situation  in  y-ray  laser  work  with  focus  on  two  areas  of  research  interest  critical  to  concepts 
for  developing  a  y-ray  laser.  Heating  effects  associated  with  upconversion  techniques  were 
discussed.  The  sources  of  inhomogeneous  broadening  which  destroy  the  Mossbauer 
Effect  were  investigated  and  techniques  available  for  restoring  the  resonance  b}  external  or 
internal  fields  were  considered  (IDA  Ptq)er  P-2083). 

In  1988  the  in-house  woik  concentrated  on  (1)  establishing  a  theoretical  model  of 
nuclear  superfluorescence  which  took  into  account  specific  characteristics  of  nuclear 
radiatitm  enussion  and  transport  that  are  important  for  superfluorescence  on  the  nuclear 
scale,  (2)  examining  in  more  detail  and  for  more  realistic  parameters  the  heating  effects 
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inherent  in  upconversion  concepts,  and  (3)  examining  the  speed  of  response  of  nuclei  to 
applied  external  fields,  as  exhibited  by  their  spectra  (IDA  Paper  P-217S). 

The  1989  effon  described  in  the  present  document  concentrated  on  developing  a 
model  of  nuclear  superfluorescence  that  included; 

1 .  Multilevel  systems, 

2.  Incoherent  pumping  at  finite  rates, 

3.  Spatial  attenuatitxi  of  the  electromagnedc  field  in  the  active  region, 

4 .  Dephasing  due  to  inhomogeneous  broadening, 

5 .  Population  depletion  due  to  competing  processes  like  internal  conversion  and 
spontaneous  transverse  emission. 

This  was  done  by  generalizing  the  Haake-Reibold  model  of  superfluorescence  based 
on  Maxwell-Bloch  equations.  The  effects  of  inhomogeneous  broadening  on  super- 
fluorescence  were  studied  using  a  theory  developed  by  J.  Eberly.  The  enhancement  of 
superfluorescence  in  inhomogeneously  broadened  systems  by  time-dependent  hyperfine 
interactions  was  investigated  using  the  theory  developed  by  M.  Blume.  We  also 
investigated  coherent  aspects  of  nuclear  emission  and  absorption  in  the  presence  of 
homogeneous  and  inhomogeneous  broadening  of  source  and  absorber. 

The  development  of  a  y-ray  laser  is  viewed  as  a  high-risk/high-payoff  undertaking. 
Ida's  involvement  focuses  on  minimizing  the  risk  and  on  striving  to  redirect  the  effort 
when  piTq>os6d  schemes  are  shown  not  to  be  feaable. 
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This  paper  generalizes  the  Haake-Reibold  model  of  supeifluorescence  to  take  into 
account  certain  effects  that,  while  negligible  for  atomic  emission,  are  important  for  the 
nuclear  case.  Comparisons  are  made  of  three  different  approaches  to  numerical 
computation  based  on  the  model,  and  a  number  of  explicit  calculations  are  made  to  illustrate 
the  behavior  of  the  emitted  pulse  under  various  conditions.  In  additirm,  the  paper  presents 
detailed  analyses  of  the  effects  of  homogeneous  and  inhomogeneous  broadening  and 
relaxation  on  superfluorescent  emission. 
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SUMMARY 


This  repent  describes  the  1989  research  effort  by  members  of  the  IDA  staff  in  the 
field  of  Y*ray  lasers.  The  work  is  part  of  a  continuing  task  in  support  of  the  Innovative 
Science  and  Technology  Office  (1ST)  of  the  Strategic  Defense  Initiative  Organization 
(SDIO).  The  development  of  a  y-ray  laser  is  a  high-risk  science  and  technology 
undertaking.  IDA  involvement  has  focused  in  large  measure  on  minimizing  the  risk  and 
attempting  to  redirect  the  program  as  quickly  as  possible  when  proposed  schemes  prove 
infeasible.  The  report  is  presented  in  three  independent  chapters. 

Chapter  I  discusses  a  generalized  version  of  the  Haake-Reibold  model  of 
superfluorescence  in  a  multilevel  system  of  emitters,  having  up  to  five  energy  levels.  The 
generalized  version  takes  into  account  finite  pumping  rates,  spatial  attenuation  of  the 
electromagnetic  field,  dephasing,  and  population  depletion  due  to  processes  such  as  internal 
conversion  and  spontaneous,  but  noncollective,  emission.  These  phenomena,  which  have 
not  been  considered  together  in  a  single  theory  heretofore,  are  important  for  the 
investigation  of  the  feasibility  of  nuclear  superfluorescence. 

Calculations  performed  with  the  Haake-Reibold  model  show  that  nuclear 
superfluorescence  may  be  less  demanding  on  the  real  nuclear  and  solid  state  properties  than 
was  previously  predicted  by  other  authors. 

Chapter  II  discusses  the  effects  of  inhomogeneous  broadening  on  reducing 
superfluotescent  emission  and  the  effects  of  time-dependent  interactions  on  enhancing  such 
emission.  The  effects  are  functions  of  time  (after  inversion)  unlike  similar  effects  in 
continuous  wave  (CW)  or  pulsed  (amplified  spontaneous  emission)  lasing. 

Chapter  m  discusses  the  effect  of  inhomogeneous  broadeiting  on  emission  and 
absorption  of  y-rays  by  nuclei.  The  phenomena  is  treated  as  incoherent  and  coherent 
processes  to  examine  claims  made  by  some  investigators  that  homogeneous  broadening 
can  be  used  to  overcome  reductions  of  resonant  absorption  produced  by  inhomogeneous 
broadening.  It  turns  out  that  in  both  cases  the  resonant  absorption  (Mdssbauer  Effect)  is 
reduced  by  a  factor  y/T,  where  y  is  the  radiative  linewidth  and  F  the  total  linewidth  of  the 
system  and  the  predicted  enhancement  does  not  occur.  The  total  linewidth  includes  the 


radiative  process,  internal  conversion,  cooperative  transitions,  and  time-dependent  fields 
that  tend  to  broaden  the  emission  line. 
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I.  INVESTIGATION  OF  NUCLEAR  SUPERFLUORESCENCE 
WITH  THE  GENERALIZED  HAAKE-REIBOLD 
STOCHASTIC  SOURCE  APPROACH 

A.  INTRODUCTION 
1 .  Background 

Superfluorescence  (SF)  and  superradiance  (SR)  are  both  examples  of  the 
cooperative,  spontaneous  emission  of  coherent  radiation  by  identical  radiators  such  as 
atoms,  molecules,  or  nuclei.  Intense,  directed  pulses  of  duration  much  shorter  than  the 
spontaneous  emission  lifetime  of  an  individual  radiator  (Ref.  1)  characterize  the  radiation, 
which  has  been  observed  experimentally  in  iix>lecular  and  atomic  systems,  e.g.,  CH3F, 
HF,  Na,  Ca,  Tl,  KCI:02,  Sr,  and  Li.  Theoretical  models  that  successfully  explain  the 
observations  have  been  developed. 

For  nuclear  lasing,  which  has  special  pumping  problems  (Ref.  2),  it  should  be 
easier  to  invert  a  population  to  an  isomeric  level  with  a  long  lifetime  than  to  one  with  a  short 
lifetime.  On  the  other  hand,  during  the  actual  lasing  process  tiie  emissitMi  lifetime  should 
be  as  short  as  possible.  Thus,  because  of  its  characteristically  fast  radiation,  some 
(Refs.  2,  3)  have  recommended  SF  as  a  possible  alternative  to  stimulated  emission  as  a 
viable  or  ptefetred  process. 

A  recent  IDA  paper  (Ref.  3)  discussed  some  aspects  of  nuclear  SF  and  compared  it 
to  lasing  and  ampUHed  spontaneous  emission  (ASE).  Results  indicate  that  electronic 
attenuation  plays  a  more  complicated  role  in  generating  SF  than  in  those  processes.  It  does 
not  lead  to  the  threshold  required  by  the  Schawlow-Townes  condition  for  lasing  and  ASE, 
but  rather  forces  a  gradual  reduction  of  emission  intensity.  The  paper's  calculations  were 
based  on  a  three-state  noodel  with  some  restrictive  conditions  that  permitted  an  effectively 
closed  form  solution. 

The  present  paper  uses  a  five-state  model  and  a  Monte  Carlo  technique  for 
numoical  computations.  The  treatment  is  more  general  and  takes  into  account  interacting 
physical  mechanisms  that  the  earlier  nuclear  SF  work  ignored. 
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2.  Comparison  of  Atomic  and  Nuclear  SF  Problems 


Nuclear  SF,  which  imposes  special  conditions  on  the  nuclei  and  their  enviraiment, 
has  not  yet  been  observed.  For  SF  certain  parameters,  comparative  values  of  which  appear 
in  Table  1,  are  unimportant  in  the  atomic  domain  but  are  critical  in  the  nuclear.  The 
parameters  in  Table  1  determine  the  gain  coefficient  K  associated  with  an  amplifying 
medium  for  typical  attxnic  and  nuclear  systems.  The  gain  coefficient  is  given  roughly  by^ 


f  r 

2jc(i  +  a)  r+r, 


(l  +  a) 


|£.n,  -Mp 


where  Or  is  the  effective  nuclear  resonance  cross-section,  and  An*  is  the  population 
inversion  achieved  by  the  pumping  process.  The  M  depends  on  the  laser  operating  riKxle 
as  discussed  by  Trammell  and  Haiuion  (Ref.  4).  For  steady  state  operation  M  =  1. 

In  Eq.  (1),  Or  is  the  maximum  nuclear  resonant  absorption  cross-section  given  by 


R  2ii(l  +  a)  (1-1- a)  ’ 

with  X  the  wavelength  of  the  photon  emitted  in  the  resonant  transition,  a  the  internal 
conversion  coefficient,  and  Pr  the  branching  ratio.  The  quantity  F  is  the  total  width  of  the 
nuclear  excited  level  decaying  with  transition  energy  Eq  ,  and  Fg  is  the  total  width  of  the 
lower  level.  Generally,  the  experimental  linewidth  is  often  much  larger  dian  F  because  of 
hoax)geneous  broadening  due  to  competing  transitions,  the  final  nuclear  state  lifetime,  and 
inhomogeneous  line  broadening.^ 

Equation  (1)  is  the  same  as  that  given  for  the  gain  coefficient  by  Baldwin  and  Feld 
(Ref.  5),  except  for  two  additions.  One,  which  takes  into  account  spatial  attenuation  of  the 
electromagnetic  field,  is  the  second  term  on  the  righthand  side.  The  other  is  the  faaor 
1/(1  -f-  a),  which  takes  into  account  inhmnogeneous  broadening. 


^  A  more  detailed  treatment  (Ref.  3)  of  SF  shows  that  the  relationship  between  the  attenuation  p  and  the 

stimulated  gain  ORAn*  is  more  complicated  than  that  depicted  in  Eq.  (1). 

The  major  contributions  to  inhomogeneous  broadening  come  from  hyperfme  interactions  such  as  the 
electric  monopole  (isommr  shift),  magnetic  dipole  (Zeeman),  and  electric  quadrupole  (electric  field 
gradient). 
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Tabitt  1.  Comparison  of  Characteristic  Atomic  (Optical)  and 
Nuclear  MOssbauer  Parameters 


Parameter 

Typical 

Optical 

Atomic 

Typical 

Nuclear 

Possible 

Modification 

Factor 

Modification 

Technique 

Energy 

^(eV) 

eV 

10^  eV 

Wavelength 

X(A) 

id*  A 

1A 

ReooHless  Fraction 

f 

1 

-  0 

(0.95) 

MOssbauer 

Effect 

Inhomogeneous 

Broadening 

a 

smal 

10^-10® 

x10~^ 

Dynamic  Reids 
Relaxation 

Competing  Transitions 
(Internal  conversion 
fornuclen 

a 

smal 

10-100 

Attenuation 

Mfc(cm~^) 

0.01 

xtO"® 

Bormtann 

Effect 

Density 

n(cm“®) 

10*°-10*®* 

10^2 

*  SF  observed  in  this  density  range  in  molecular  and  atomic  systems  (Ref.6). 


The  inhomogeneous  broadening  parameter  a  can  be  10^  or  higher  for  istmieric 
nuclei  of  interest  to  the  f-ray  laser  concept  under  consideration.  The  parameter  g^  is  the 
coupling  fiikctor  for  the  nuclei  to  the  electromagnetic  field,  which  for  amorphous  samples  is 
equal  to  1,  and  the  parameter  p  is  the  linear  extinction  coefficient  for  nonresonant 
absoqjtkm  or  inelastic  scattering  of  jdiotons  primarily  1^  the  photoelectric  effect 

Hrst,  it  is  clear  from  Eq.  (1)  diat  the  nuclei  participating  in  lasing  have  to  belong  to 
a  Mdssbauer  isotc^  (f  >  0),  as  even  simple  nuclear  resonant  scattering  is  not  possible 
without  this  condition  being  fulfilled.  Preferably,  the  Mbssbauer  isotope  should  be  in  a 
medium  such  that  one  has  a  high  recoilless  fractitm  (f  «  1.0).  Thus,  the  pumping 
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mechanism  must  be  capable  of  achieving  the  required  nuclear  inversion  (An*)  without 
destroying  the  Mbssbauer  Effect  in  the  crystal.  For  direct  pumping,  the  pumping  power 
requirement  for  inverting  the  nuclear  population  of  a  given  isotope  is  inversely  proportional 
to  the  lifetime  t  of  the  nuclear  state.  Thus  in  Table  1  we  see  the  recoil  (f  ■»  0) 
inhomogeneous  broadening  (a  >  lO^),  internal  conversion  [high  (a  >  10)  for  low  energy 
Y-emission]  and  high  attenuation  (p.  >  10)  combine  to  reduce  the  gain  coefficient  and  the 
feasibility  of  lasing  and  SF  for  nuclear  systems,  whereas  the  possibility  of  high  densities 
(n  »  10^2)  enhances  the  gain  and  feasibility  of  lasing  and  SF  for  nuclear  systems. 

3.  SF  Pulse  Characteristics 

Generally  speaking,  the  emission  of  an  SF  pulse  requires  the  preparation  of  an 
inverted  population  of  identical  radiators.  Experimentally,  the  emitted  pulses  are 
characterized  by  a  pulse  width  xgR  and  a  delay  time  tD>  following  the  inversion  of  the 
populatimi.  From  simple  theory,  the  relationship  between  xgR  ‘^D  &nd  the  density  of 
cooperating  radiators  p,  tire  natural  radiative  lifetime  xo,  the  wavelength  of  tire  emitted  light, 
X,  and  the  cavity  length  /  can  be  obtained.  In  terms  of  these  parameters 


Xp  =  -i-x^lnN, 

and  the  pulse  intensity  Is(t)  is  given  the  "inverted  pendulum  expression” 

Is(t)  =  1  N/x^  sech^['^  ■  V] 

(Ref.  7),  where  t  is  measured  from  the  time  of  the  (instantaneous)  inversion. 

Qearly,  to  describe  experimental  conditions  more  realistically  it  is  necessary  to 
invttite  trxxe  detailed  theories,  which  involve  the  parameters  in  Table  1  and  Eq.  (1). 

4.  Historical  Approaches  to  Nuclear  SF  and  Objectives  of  this  Paper 

This  paper  discusses  an  SF  model  that  takes  into  account;  (1)  competing 
transitions,  such  as  internal  conversion;  (2)  transport  phenomena  causing  electromagnetic 
spatial  attenuation;  (3)  finite  pumping  times  for  the  preparation  of  inversion  by  incoherent 
sources;  (4)  effects  of  inhtmxrgeneous  broadening  on  pulse  etrussion.  These  phenomena 
lead  to  conditions  for  SF  emission.  Although  characterized  by  distinct  parameters,  they 


must  be  treated  together  in  the  same  theory  since  the  interplay  between  the  different  effects 
that  they  produce  is  too  intricate  to  permit  than  to  be  ueated  separately. 

Some  effects  have  been  considered  separately  by  other  authors,  however. 
TrammeU  and  Hanix>n  (Ref.  4)  addressed  the  feasil^ty  of  nuclear  SF  and  compared  SF  to 
nuclear  "pulsed"  lasing  or  amplified  spontaneous  emission  (ASE).  They  considered  the 
effects  on  SF  of  inhomogeneous  broadening  and  attenuation  separately  and  introduced 
these  parameters  in  an  od  hoc  fashion  in  the  discussion.  Their  conclusions  differ 
substantially  from  ours,  in  two  important  features.  In  our  results,  the  effect  of  the 
attenuation  is  more  gradual,  and  does  not  lead  to  a  threshold  effect  as  required  by  the 
Schawlow-Townes  condition  (see  Ref.  3)  for  lasing.  Also  in  our  results  the 
inhcnnogeneous  broadening  leading  to  a  dephasing  effect  is  time  dependent  The  effect  is 
small  at  eariy  times  after  inversion  and  increases  with  time.^  Both  of  these  effects  when 
treated  together  in  our  theory  produce  a  less  restrictive  condition  on  the  feasibility  of 
nuclear  SF.  Baldwin  and  Feld  in  a  more  recent  paper  (Ref.  S)  using  semiclassical  theory, 
discuss  SF  emission,  and  derive  minimum  requirements  for  an  e]q>erimental  demonstration 
of  phenomena.  They  do  not  consider  inhomogeneous  broadening,  and  grossly 
overestimate  die  effect  of  attenuation. 

We  plan  to  use  the  more  complete  thetny  developed  in  this  paper  to  reevaluate  the 
feasibility  of  demonstrating  experimentally,  with  known  isomers  and  available  pumping 
sources,  the  occurrence  of  nuclear  SF.  The  model  described  here  assumes  that  a  single, 
endfire  mode  exists  in  an  active,  acicular-shaped  medium  (i.e.,  a  medium  shaped  like  a 
long  cylinder).  Geometrical  considerations  leading  to  the  treatment  of  more  general  shapes 
(other  than  acicular)  will  be  left  for  future  discusnons. 

5.  SF  Theory  for  Atomic  Systems 

In  the  introduction  to  Ref.  7  the  authors,  Bonifacio  and  Lugiato,  distinguish 
between  SR  and  SF.  Both  terms  refer  to  the  phenomenon  of  cooperative  radiation  by  a 
system  of  identical  emitters,  which  in  their  case  are  two-level  atmns. 

However,  SR  occurs  when  a  odierent  pulse  excites  die  system  into  a  collective  state 
with  a  macroscopic  electric  dipole;  the  supenadiant  pulse  is  then  attributed  to  emission  by 
this  dipole.  On  the  other  hand,  SF  occurs  when,  after  some  normal  fluorescent  emission  to 


The  rate  is  die  same  but  there  is  an  ciqioneiuial  decrease  in  effect 
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start  the  process,  with  no  macroscopic  dipole  present  initially,  uncorrelated  excited  emitters 
spontaneously  become  correlated  and  then  radiate  cooperatively.  At  no  time  is  there  a 
macroscopic  dipole  present. 

For  the  Bonifacio  and  Lugiato  thecxetical  model  of  SF,  the  active  volume  containing 
the  emitters  is  limited  to  an  acicular  shape:  many  wavelengths  long  and  with  a  cross- 
section  narrow  enough  to  make  the  Fresnel  number  associated  with  diffraction  of  the 
radiation  at  one  end  neariy  equal  to  1.  This  condition  limits  the  electromagnetic  field  within 
the  volume  to  a  pair  of  modes,  characterized  to  a  good  approximation  by  left  and  right 
nooving  plane  waves.  The  authors  assume  that  these  modes  are  completely  independent 
and  noninteracting,  thereby  reducing  the  analysis  inoblem  to  one  for  a  single  nxxle. 

In  their  analysis  Bonifacio  and  Lugiato  consider  both  Markovian  and  non- 
Markovian  cases  of  the  interaction  between  the  local  electromagnetic  field  and  the  atoms. 
For  the  non-Markovian  case,  in  which  a  certain  number  of  the  emitted  photons  remain 
within  the  active  volume  long  enough  to  reexcite  some  depleted  atoms,  the  model  predicts 
an  experimentally  observed  ringing  in  the  emitted  pulse.  For  tiie  Markovian  case  the  nKxiel 
predicts  that  ringing  will  not  occur  because  all  emitted  photons  leave  the  volume  before  they 
can  interact  with  any  of  the  atoms. 

The  Bonifacio  and  Lugiato  model  is  based  entirely  on  the  use  of  a  density  matrix 
that  is  a  function  of  time  but  has  no  spatial  dependence.  Therefore,  it  does  not  take 
propagation  into  account  directly.  Nevertheless,  it  is  able  to  predict  ringing,  whidi  to  some 
extent  is  due  to  prt^agation. 

In  this  context  it  has  been  referred  to  as  an  average  field  model.  Some 
(cf.  Refs.  8,  9)  contend  that  an  average  field  model  is  incapable  of  predicting  such 
phentmiena  as  ringing  correctly,  arguing  that  it  is  not  only  necessary  to  treat  propagation 
explicitly  but  also  to  include  certain  three-dimensional  effects  as  well  (e.g.,  field  variations 
in  the  lateral  directions). 

Theoretical  models  that  do  treat  propagation  three  dimensionally  are  either 
inappropriate,  if  the  medium  is  crystalline,  because  they  assume  that  the  etrutters 
are  distributed  randomly  over  some  region  (cf.  Ref.  10);  or  else  they  are  semiclassical 
(cf.  Ref.  1)  and  therefore  do  not  deal  completely  with  the  iititial  stage  of  superfluorescence. 
On  the  other  hand,  the  medium  considered  specifically  by  Bonifacio  and  Lugiato  is,  in  fact, 
a  crystal,  and  the  Bonifacio  and  Lugiato  model  is  strictly  quantum  mechanical  in  nature. 
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If  the  emitters  are  nuclear  rather  than  atomic,  which  is  the  case  of  interest  here,  the 
multilevel,  <Hie-dimensional  SF  model  of  Haake  and  Reibold  (Ref.  11)  may  be  more 
suitable  than  any  of  these.  The  effect  of  a  crystalline  medium  should  be  particularly 
important  when  the  emitters  are  nuclei  because  the  spacing  between  them  will  generally  be 
larger  than  a  lesoiuint  wavelength,  whereas  when  they  are  atoms  the  opposite  is  true.  The 
Haake  and  Reibold  model  covers  at  least  one-dimensional  propagation  completely,  and, 
since  quantum  electro-dynamics  plays  a  majOT  role  in  its  derivation,  it  should  also  cover  the 
initial  stage  of  SF  adequately.  The  generalized  version  of  the  model  defined  here  takes  into 
account  certain  effects,  such  as  losses  due  to  internal  conversion,  that  are  negligible  in  the 
atomic  case  but  are  in:qx>itant  in  die  nuclear  case. 

B .  HAAKE-REIBOLD  STOCHASTIC  SOURCE  APPROACH 

The  Haake-Reibold  theoretical  model  of  SF  described  in  this  section  is  a 
generalization  of  the  one  introduced  by  those  authors  in  Ref.  11.  For  its  derivation 
Ref.  1 1,  which  treats  multilevel  systems,  refers  to  earlier  papers  (e.g..  Ref.  12^)  that  deal 
with  two-state  systems  using  the  same  approach.  These  references  give  more  or  less 
detailed  derivations  of  the  two-state  version  only;  however,  as  Haake  and  Reibold  remark, 
diey  lay  out  all  concepts  needed  for  generalizing  to  the  multilevel  case. 

Both  versions  use  the  Heisenberg  picture  and  both  approximate  the  discrete  set  of 
emitters  as  a  cmitinuous  distribution  in  space.  In  effect,  the  polarization  and  occupation 
numbo*  (invasion)  operators  associated  with  the  different  energy  levels  become  fields 
represented  by  continuous  functions  of  distance  and  time.  The  Heisenberg  equations  of 
nootion  are  formally  similar  to  the  Maxwell-Bloch  equations  of  semiclassical  noodels  such 
as  that  presented  in  Ref.  1. 

Howeva,  in  deriving  die  semiclassical  equations  (using  a  classical  electromagnetic 
field  and,  to  take  the  emitters  into  account  quantum  mechanically,  a  nonrelativistic. 


*  Reference  12  proves  that  a  classically  stochastic  noise  source  can  rqilace  the  initial  quantum 
electrodynan^  vacuum  fluctuations  that  are  actually  reqxmsible  for  starting  the  superfluorescmice 
process.  This  makes  it  possible  to  formulate  an  equivalent  semiclassical  initial  value  problem,  the 
acriution  of  which  provides  the  temporal  history  of  the  radiated  electromagnetic  pulse.  The  proof  goes 
even  firrther,  demonstratmg  that  the  noise  source  is  Gaussian  distributed  and  white. 

Reference  12  also  derives  a  one-dimensional  SF  model  characterized  by  MaxweO-Bloch  equations  in 
which  die  polarization  and  inversion  operators  are  continuous  functions  of  space  as  well  as  time.  The 
derivation  predicts  homogeneous  broadening,  giving  an  explicit  value  for  the  natural  lifetime  of  the 
S]fStem  of  emitters.  That  is,  the  cotiesponding  nature  decay  rate  is  not  an  ad  hoc  assumption  in  this 
model, 
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Schroedinger  wave  equation).  Ref.  1  at  first  obtains  the  equations  in  terms  of  expectation 
values  of  the  quantum  mechanical  observables:  the  induced  polarization  and  the  inversion. 
These  quantities  are  continuous  functions  of  time  at  each  emitter  and  are  coupled  between 
the  different  emitter  locations  only  through  the  electric  field.  Reference  1  replaces  them 
with  polarization  and  inversion  densities,  which,  in  effect,  are  classical  fields  also  coupled 
to  the  electric  field  but  depending  continuously  on  both  space  and  time. 

After  an  inversion  of  a  sufficiently  large  number  of  emitters  to  an  unstable  energy 
level,  quantum  fluctuations  of  the  itutial  vacuum  state  electromagnetic  field  are  generally 
regarded  as  responsible  for  triggering  SF.  In  order  to  avoid  the  complications  of  a 
quantum  field  theoretical  analysis.  Ref.  1  introduces  an  ad  hoc  polarization  noise  source  to 
start  the  process.  Reference  12  proves  that  a  classically  stochastic  polarization  source, 
properly  chosen,  can  in  fact  substitute  for  the  quantum  fluctuations,  thereby  justifying  in 
part  the  less  rigorous  approach  of  Ref.  1. 

The  stochastic  elements  of  the  source  assumed  in  Ref.  1  are  terms  with  random 
phases  that  when  summed  have  little  effect  on  the  early  electric  field  intensity,  which  can  be 
estimated  by  neglecting  iK>nlinear  interacticms.  The  nonlinear  interactitms  are  small  because 
tiiey  are  quadratic  in  quantities  that  vaiush  in  the  limit  as  time  approaches  zero. 

Reference  1  uses  the  results  of  such  a  calculation  for  a  time  near  zero,  but  after  the 
noise  source  becomes  negligible,  to  obtain  a  value  of  the  electric  field.  The  quantity  thus 
acquired  is  then  available  as  an  iiutial  value  fen*  an  irutial  value  problem  associated  with  the 
Maxwell-Bloch  equations,  the  solution  of  which  will  determine  the  field  at  later  times  when 
die  nonlinear  interactions  ate  no  longer  negligible. 

The  Haake  and  Reibold  models  also  introduce  a  classically  stochastic,  Gaussian 
white  noise  source  to  replace  the  quantum  field  fluctuations  as  a  means  of  starting  the 
collective  spontaneous  emisssion  that  characterizes  superfluorescence,  but  use  the  Ref.  12 
result  for  this  purpose.  In  addition,  the  introduction  of  other  source  terms  of  the  same  type 
provides  the  multilevel  model  with  incoherent  depletion  effects  at  lower  energy  levels. 

Reference  11  determines  the  second  order  statistics  of  the  electric  field  and 
polarization  at  early  times  when  the  Maxwell-Bloch  equations  are  effectively  linear.  Then, 
treating  the  Maxwell-Bloch  equations  as  classical,  it  uses  the  results  to  generate  a  statistical 
distribution  of  initial  values  for  calculating  a  cesresponding  distribution  of  solutions  that  ate 
valid  in  the  nonlinear  regime.  Essentially,  this  is  a  Monte  Carlo  approach  to  estimating 
various  features  of  the  emitted  radiation,  such  as  the  delay  time  of  the  radiated  pulse. 


Although  restricted  to  one  spatial  dimension,  the  Haake  and  Reibold  models 
account  explicitly  for  propagation  of  the  emitted  electromagnetic  field  as  long  as  the 
geometrical  configuration  of  the  emitters  limits  the  radiation  to  an  endfire  mode.  For  this 
purpose  the  assumed  shape  of  the  active  region  is  acicular  with  an  associated  Fresnel 
number  of  the  order  of  one. 

The  effect  is  to  restrict  the  field  to  plane  waves  propagating  to  the  left  or  to  the  right 
in  the  Itmgitudinal  direction.  The  models  assume  that  these  waves  do  not  interact,  i.e.,  that 
it  is  (Mily  necessary  to  consider  one  of  them  in  calculating  the  radiation  at  a  single  end  of  the 
active  region. 

The  most  general  multistate  system  of  interest  for  analyzing  the  nuclear  SF  is  one 
with  five  states  for  which  the  energy  diagram  is  depicted  in  Fig.  1.  Pumping  is  modeled  by 
the  transition  from  the  storage  level  4  to  this  upper  lasing  level  3.  We  allow  for  two  lasing 
transititms  (3  to  1  and  3  to  2)  and  the  p(»sibility  of  interference  between  them.  Level  2  can 
depc^ulate  to  level  1  and  level  1  to  level  0. 


0 


Figure  1.  Energy  Level  Structure  Assumed  in  the  Calculations  Described  In  the 
Paper.  The  two  competing  SF  transitions  are  between  upper  state  3  and  lower 
states  1  and  2  and  are  Indicated  by  gi  and  g2.  respectively.  Level  4  populates 
level  3,  providing  the  Inversion.  The  decay  rates  y,  Ta,  r,,  and  r-i  are  given 

by  *  ‘^92}  '<21  respectively.  Other  depopulation  transitions  are 

generally  allowed,  but  we  Indicate  only  those  of  Interest  at  presem. 

The  corresponding  Haake  and  Reibold  Maxwell-Bloch  equations  are 

=  (2.,) 

^*^3  “”(®3I^l'*’®32*S2'*’^l  ^2  **32)” ^3  ^3  ~ 5.  (2-b) 


1-9 


N,  =  -^  (e>I,  ^  E-,  r;,  ) + r^  -  r,  N,  +  ( 1  -  p)  r,  Nj + 


1. 

dt 


—  N  =r  N 
dt  ^0  ^  ri 


dt  *^32  “  (^3  ^2)^  ^1  *^21  2  [  ^3  ■‘'  ^2"^  ^4.2]*^32  '*'  ^32 

I<'(N3-N0^1  *^2,  4[^3^ 

di  *^21  “  ^31  *^32  '‘'  ^32  ^21”  1  [  ^1  '*'  ^2  ]*^21 '‘’  ^21 

lr4=8K4^^4<'='*2>  • 


(2-0 

(2-d) 

(2-0 

(2-f) 

(2-g) 

(2-h) 

(2-i) 


which  are  a  generalized  version  of  the  analogous  equations  introduced  in  Ref.  1 1  in  which 
the  population  number  operators  are  Ni,  the  photon  field  operators  are  Eij^,  the  nuclear 
polarization  operators  are  R^ij  and  the  natural  decay  rates  are  Pi.  Note  that  Pi,  associated 
with  a  transition  from  the  i^  energy  level,  is  a  decay  rate  that  is  not  due  to  the  collective 
emission  process;  rather,  Pi  takes  into  account  decay  due  to  other  spontaneous,  but 
noncollective,  processes,  such  as  internal  conversion.  Appearing  in  Eqs.  (2-c)  through 
(2-i)  are  four  additional  parameters:  p.,  which  is  a  spatial  attenuation  coefficient,  and 
P^i  (t=lt  2)  which  are  homogeneous  line-tooadening  decay  constants  and  act  as  dephasing 
parameters,  and  P  which  is  a  branching  ratio  for  emissions  from  state  3  to  states  1  and  2. 
In  addition,  the  time  dependence  of  the  coupling  parameter  g,  which  in  the  Haake-Reibold 
theory  is  constant,  introduces  the  inhomogeneous  broadening  effect  on  the  SF;  it  is 
assumed  to  be  exponential:  gi'(t)  =  gi  cxp[-l/2  Peit],  where  Pei,  which  is  the 
inhomogeneous  broadening  decay  constant,  is  also  the  decay  rate  for  the  coupling 
parameter.  Both  P^i  and  Pei  can  be  written  as  proportional  to  the  natural  line  widths 
broadening  parameters  b  and  a  (assumed  for  simplicity  to  be  the  same  for  both  transitions 
3  2  and  3  -♦  1),  so  that  P^i  =  bP/  and  Pei  =  aPi^  where  Pi^  =  Pi  P3.  The  effect  on  SF 

of  the  inhomogeneous  broadening  through  the  ^me  dependence  of  the  coupling  parameter  g 
has  been  studied  by  several  investigators  (Refs.  7, 13).  Chapter  n  of  this  report  shows 
how  the  effect  of  inhonoogeneous  broadening  can  be  reduced  through  relaxation  or  time- 
dependent  hyperfine  interactions. 

The  inhomogeneous  terms  of  the  form  ^  in  Eqs.  (2)  are  noise  sources  that  cause 
incoherent  depletion  at  the  various  energy  levels.  The  source  ^  provides  for  incoherent 
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pumping  of  energy  level  4,  from  which  energy  decays  at  the  same  rate  to  level  3,  causing 
the  population  inversion  that  eventually  results  in  SF.  The  inhomogeneous  terms  of  the 
form  represent  the  noise  sources  that  the  model  substitutes  for  the  initially  important 
quantum  field  fluctuations  pumping  and  spontaneous  decay  effects  on  the  SF  triggering 
process.^ 

Reference  1 1  assumes  that  the  fields  ^  and  ^32,  although  quantum  mechanical 
operators,  are  stochastic  in  the  classical  sense  with  the  following  statistical  properties. 
They  are  statistically  independent,  white,  and  Gaussian  with  zero  means  and  the  second 
moments 

<  ^32'"  ^32“^  >  =  <  432”  ^32"  >  =  <  ^32~  ^32'^  >  =  0.  (3) 

<  ^  (x,  t)  ^4  (x’,  t')  >  =  <  432"^  (x,  t)  432"  (x',t')  > 

=  (1/N)  <  N4  (t)  >  y5(x  -  X’)  5  (t  -  f)  (4) 

=  (l/N)e-1^76(x-x’)5(t-t’)  . 

A  comparistxi  of  Eq.  (3)  with  Eq.  (4)  verifies  that  the  quantities  432'‘'  and  432“ 
in  fact,  operators  rather  than  c-numbers,  since  they  do  not  commute.  The  term  y  in 
Eq.  (4)  is  dimensionless  and  is  defined  by  the  ratio  of  the  superfluorescence  time  to  the 
natural  decay  time  from  level  4  to  level  3: 

Y=tsF^T43  .  (5) 

The  quantity  N  is  the  number  of  radiators. 

Although  the  polarization  sources  432  and  the  population  source  44  are  statistically 
independent,  as  indicated  in  Eqs.  (3)  and  (4),  Ref.  1 1  assumes  that  the  parameters  defining 


^  The  stochastic  n(^  operatots  are  defined  by  their  statistics.  Polder,  Shuurmans,  and  Vrehens  (Ref.  12) 
derived  the  classical  noise  equivalent  sources  from  a  quantum  mecimical  model  They  identified  three 
distinct  components  of  the  EM  field,  namely,  (1)  vacuum  fluctuations,  Evac>  (2)  E^a  os  dipole-dipole 
intoaction  b^een  difierent  atoms  and  (3)  the  self  force  or  radiative  reaction,  E„^.  Initially,  the  only 
nonzero  component  is  Evac-  fhis  sets  up  correlations  in  the  polarizations  of  atoms  which  produces  a 
nonzero  Edd.  The  Evac  wave  vectors  conreqxMiding  to  the  appropriate  geometry  of  the  active  region 
(cavity)  set  up  the  polarization  appropriate  for  cooperative  emission  with  those  wave  vecuns.  This 
results  in  the  observed  SF  emission.  The  E,t  field  gives  rise  to  the  natural  width  of  the  emission  line 
only.  [This  interpretation  arises  naturally  from  the  anti-normal  ordering  chosen  by  Polder  et  al.  For  a 
discussion  of  the  influence  of  ordering  of  operators  on  interpretive  questions,  see  P.W.  Milonni  and 
W.A.  Smith,  Pfiys.  Rev.  A.  (1975).]  Haake  and  Reibold  (Ref.  11)  use  this  theory  to  model  noise 
sources  which  trigger  SF.  Their  noise  source  is  due  to  the  incoherent  pumping  mechanism  with  time 
constant  y  or  the  natural  qxmtaneous  decay  from  the  excited  state  with  time  constant  r3.  In  the 
present  development  we  use  the  noise  source  from  the  incoherent  pumping  mechanism,  as  described  in 
Eqs.  (4),  for  the  SF  triggering. 
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their  probability  distributions  (ix.,  the  first  and  second  moments)  are  identical.  This  is 
presumably  justified  by  the  fact  that  the  same  radiation  field  reservoir,  consisting  initially  of 
the  vacuum  state  electric  field,  generates  both.  The  form  of  the  source  terms  on  the  right- 
hand  side  of  Eq.  (2)  is  consistent  with  the  stochastic  source  derived  in  Ref.  12. 


C.  COMPUTER  CALCULATIONS,  PART  A 

The  computer  calculations  to  be  described  here  were  done  for  only  one  SF 
transition  (3  1),  the  other  transition  was  supressed.  The  calculations  served  two 

purposes:  to  check  the  codes,  and  to  leam  about  the  physics  of  nuclear  SF  predicted  by  the 
theory  of  Chapter  n.  We  first  wanted  to  check  our  codes  by  comparing  our  results  with 
those  obtained  by  other  investigators  and  also  to  relate  the  calculated  pulse  shapes  to 
physically  meaiungful  parameters.  We  (1)  investigated  the  dependence  of  the  pulse  shapes 
on  the  statistics  and  time  dependence  of  the  source  terms  ^i;  (Z)  compared  the  full  Monte 
Carlo  calculations  with  the  approximate  result  fix>m  the  closed  form  solution  obtained  with 
the  "triple  integral"  approach  already  introduced  and  described  in  an  earlier  report  (Ref.  3); 
(3)  compared  their  results  with  the  average  peak  calculation  described  by  Haake  and 
Reibold  (Ref.  11).  We  also  studied  the  variation  of  the  individual  pulses  obtained  with 
different  seeds  for  the  random  number  generator  used  to  provide  the  noise  source 
characteristics  that  simulate  die  quantum  fluctuatitms.  All  of  these  results  are  described  in 
this  section. 

In  Section  IV  we  describe  the  results  obtained  with  the  Monte  Carlo  calculations 
designed  to  investigate  and  predict  the  SF  pulse  dependence  on  experimentally  achievable 
parameters  that  could  be  used  to  select  isomer  candidates  for  nuclear  SF  experiments.  The 
parameters  studied  were  (1)  the  linear  attenuation  coefficient,  p(cm~l),  (2)  the  lasing  level 
spontaneous  lifetime  F  (s~^),  including  both  radiative  emission  and  internal  conversion, 
(3)  the  inhomogeneous  broadening  parameter  p  giving  tire  effective  increase  in  linewidth  of 
the  total  system  in  units  of  the  natural  linewidth  and  the  effective  number  of  cooperating 
nuclei  N. 


All  of  the  calculated  results  are  presented  in  normalized  units  consistent  with  the 
usage  in  Eqs.  (2-a)  through  (2-i).  In  these  units,  time  is  given  in  units  of  the  SF  time  tsf 
defined  as 


3pX^l 
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where  to  is  the  natural  lifetime  of  a  solitary  nucleus,  p  is  the  density  in  the  active  volume, 
X  is  the  wavelength  of  the  emitted  radiaticm,  and  the  length  is  normalized  to  1,  the  length  of 
the  cylindrical  active  region.  In  Fig.  2  we  show  a  typical  SF  pulse  and  define  some 
parameters.  The  phenomenon  giving  rise  to  multiple  pulses  (ringing)  has  been  discussed  in 
the  literature  (Refs.  3, 7).  We  will  be  mostly  concerned  with  single  pulses.  In  Fig.  3  we 
ccm^are  the  results  obtained  with; 

(1)  the  Monte  Carlo  solution  of  Eqs.  (2-a)  through  (2-i), 

(2)  the  "triple  integral"  solution  described  in  Ref.  12, 

(3)  the  average  peak  result  obtained  by  Haake  and  Reibold  (Ref.  1 1)  and  described 
in  Appendix  B. 


Figure  2.  Characterittics  of  SF  Pulaaa.  Ganarally,  one  obsarvea  one  or  more 
pulses.  For  a  single  pulse  (or  the  first  pulse)  Indicated  by  the  heavy  continuous 
line,  we  define  a  delay  time  xoi  and  a  width  xw  In  units  of  xsf  >nd  an  Intensity  lo. 
In  the  Inverted  pendulum  model  (Ref.  7),  xw  ^  xsf>  K  Is  observed 
experimentally  In  general  that  xw  >  tsF  nnd  this  Is  also  obtained 
from  the  more  sophisticated  models.  The  other  pulses,  Indicated 
by  the  dashed  line.  If  present,  have  their  own  characteristics. 
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N  (number  of  radiating  systems) 


Figure  3.  SF  Delay  Time  to  versus  Number  of  Cooperating  Systems  N. 


The  triple  integral  solution  thus  obtained  is  based  on  an  averaging  procedure  and 
average  field  approximation.  Its  main  benefit  is  that  the  key  parameters  such  as  N,  F,  and 
p.  are  explicitly  exhibited  in  the  final  analytical  result  and  thus  their  effect  on  the  evolution 
of  the  pulse  shape  can  be  easily  followed  and  understood  in  physically  meaningful  terms. 
Averaging  over  the  noise  sources,  using  the  assumed  statistics  and  definitions  of  the 
correlations  as  they  appear  in  the  analyses,  also  obviates  the  need  for  Monte  Carlo 
numerical  techniques.  The  disadvantage  is  the  loss  of  generality.  The  average  peak  delay 
time  was  calculated  from  an  expression  (procedure)  derived  in  Ref.  11. 

The  three  curves  giving  the  delay  time  to  as  a  fiinctitm  of  the  number  of  cooperating 
nuclei,  N,  for  the  three  different  calculations  differ  about  40  percent  at  the  low  values  of 
N  and  less  than  25  percent  at  the  high  values  of  N  in  Fig.  3.  This  is  satisfactwy  agreement 
considering  the  different  assumptions  made  in  arriving  at  the  results.  The  Monte  (Zarlo 


result  was  obtained  1^  assuming  a  standard  deviation  varying  with  time  as 


in 


the  Gaussian  noise  distribution  consistent  with  Ref.  10.  This  condition  was  imposed  on 
the  randtxn  number  generator  statistics  as  described  in  Appendix  A. 


In  the  triple  integral  (TI)  result  the  noise  terms  were  rq)laced  by  their  correlations 
after  averaging  over  space.  One  would  expect  a  shorter  delay  because  of  the  loss  of  spatial 
variabiliQr  and  reduction  of  the  effect  of  the  lower  fields  at  the  left  end  of  die  active  region. 
This  can  be  seen  more  clearly  whoi  we  look  at  the  spatial  distribution  of  fields  in  the  active 
volume.  The  difference  in  results  obtained  with  the  TI  calculation  and  the  average  peak 
calculation  is,  not  surprisingly,  small,  about  5  percent  at  the  high  values  of  N  and 
negligible  at  the  low  values. 


In  Fig.  4  we  show  the  variation  of  the  Monte  Cario  results  for  different  values  of 
the  variance  assumed  in  the  noise  distribution.  Higher  variance  means  nx>re  of  the  larger 
noise  components  in  the  fluctuations,  consequently  shorter  delay  times.  The  average  peak 
and  TI  results  of  Fig.  3  fall  somewhere  around  the  value  of  variance  in  the  noise  of  10. 
The  noise  source  in  all  of  the  Monte  Carlo  calculations  varied  exponentially  with  time  with 
the  exception  of  the  results  shown  in  Fig.  S.  For  these  calculations  in  addition  to  the 
exponential  variation,  a  square  pulse  and  a  step  function  pulse  were  used  as  shown  in  the 
insert  All  of  the  pulses  had  the  strength  of  the  exponential  at  t »  0  with  the  square  pulse 
mxizero  between  t  =  0  and  t2,  the  second  time  stq>  in  the  calculation.  The  step  pulse  gives 
the  same  results  as  the  exptmental  for  low  values  of  the  pumping  rate  y,  whereas  die  square 
pulse  gives  a  longer  delay  at  these  pumping  rates.  Note  that  the  SF  pulse  delay  Xd> 
indicating  the  appearance  of  the  SF  pulse,  is  on  the  (xder  of  the  punqiing  time  (Y~*).  Thus 
the  noise  source  has  a  longer  time  to  act  on  and  affect  the  radiators.  For  larger  pumping 
rates  the  exponential  and  square  noise  pulses  give  identical  results  and  the  produced  by 
the  step  pulse  also  approaches  the  same  result 


Finally,  in  this  basic  computer  study  we  looked  at  the  variation  in  the  individual 
pulses  as  the  seed  of  the  random  number  generator  was  varied.  Figure  6a  shows  the 
individual  pulses  and  Fig.  6b  some  average  results.  For  our  work,  averaging  25-100 
pulses  seemed  sufficient  The  statistics  of  the  variation  in  the  pulses  as  represented  by  td, 
Xw  >  and  lo,  ate  discussed  by  Haake  and  Reibold  (Ref.  1 1). 


1-15 


Pulse  Time 

(In  units  of  tsF) 


Graph  llc/MAC 


Figure  4.  Delay  Time  tq  versus  Number  of  Cooperating  Systems  N.  The 
numbers  on  the  curves  represent  different  values  of  the  noise  variance. 


Figure  5.  SF  Pulse  Emission  as  a  Function  of  the  Time  Variations  in  the 
Strength  of  the  Noise.  Square  puiae,  exponential  pulse,  and  step  pulse 
time  variations  In  the  noise  strength  were  used  as  shown  In  the  Inset. 


D.  COMPUTER  CALCULATIONS,  PART  B 


In  this  chapter  we  first  discuss  the  pulse  characteristics  as  a  fimction  of  the  pumping 
rate  y,  the  linear  attenuation  coefficient  |X,  and  the  inhomogeneous  and  homogeneous 
broadening  parameters  a  and  b  for  different  numbers  of  cooperating  nuclei.  We  then 
present  some  results  which  show  both  the  spatial  distribution  of  fields  in  the  active  region 
and  the  temporal  pulse  evolution  under  different  conditions.  The  spatial  distributions  can 
be  displayed  on  the  terminal  screen  of  a  cmnputer  as  a  series  of  time-delayed  curves 
simulating  a  movie  of  the  phenomena  occurring  in  the  active  region. 

In  Fig.  7  we  show  the  pulse  emission  delay  time  as  a  function  of  the  pumping  rate 
for  three  values  of  the  cooperation  number  (N  =  1(P,10®,  10^^).  In  units  of  Tsf  for  larger 
N  the  delay  time  increases,  but  in  physical  units  (seconds)  the  delay  time  actually  decreases 

because  tsp  goes  like  .  If  the  pumping  time  is  essentially  over  before  the  pulse  starts, 
the  pumping  rate  has  no  effect  on  the  pulse  shape. 


Pumping  Raleyllgp 


Figure  7.  SF  Pulse  Emission  Delay  Time  to  as  a  Function  of  the  Pumping  Rate  y 
for  Differem  Vaiues  of  the  Cooperation  Number  N.  Note  that  aithough  the 
deiay  time  is  greater  for  iarger  N  when  expressed  in  units  of  tsf  if  ia  actuaiiy 

smaiier  when  expressed  in  seconds  because  tsp  goes  iike  ^ . 


We  have  already  discussed  die  nondueshold-like  effect  of  the  attenuation  coefficient 
^  on  the  SF  phemmiena  (Ref.  3).  The  attenuation  shifts  the  peak  to  longer  times.  If  there 
are  no  dephasing  mechanisms  reducing  the  dipole  being  built  up  by  the  quantum 
fluctuations  (noise  term  the  SF  peak  will  decrease  but  not  completely  disappear.  This 
will  htqipen  unless  the  inversion  is  destroyed  by  some  depopulation  process,  like  internal 
conversion  or  some  other  competing  transition.  This  gradual  decrease  in  the  peak  due  to  p 
is  characteristic  of  SF.  A  stimulated  emission  process,  on  the  other  hand,  shows  a 
threshold  beyond  which  no  gain  occurs. 

Figures  8  and  9  exhibit  these  trends,  with  Fig.  8  showing  the  increase  in  the  delay 
time  as  p  increases,  and  Fig.  9  showing  the  decrease  of  the  peak  intensity,  which  will 
eventually  reach  the  value  of  uncooperadve  sptmtaneous  emission  as  p  is  increased. 


UnMT  AtlMNiallon CoaMclMil  |i  (on' ) 

Figur*  8.  SF  PuIm  Emission  Dslsy  Tims  to  ••  a  Function  of  the  Linear 
Attenuation  Coefficlem  p  for  Different  Values  of  the  Cooperation  Number  N. 
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Figurt  9.  SF  PulM  Emission  Intsnslty  lo  as  s  Function  of  ths  Linear  Attenuation 
Coefficient  \i  for  Different  Values  of  the  Cooperation  Number  N. 


We  now  turn  to  the  line-broadening  mechanisms.  Homogeneous  broadening 
appears  in  the  Maxwell-Bloch  equations  only  through  the  dephasing  rate  =  bFs,  where 
b  is  the  homogeneous  broadening  parameter  and  Fs  the  natural  decay  rate  of  level  3. 
Inhomogeneous  broadening  is  taken  account  of  in  the  Maxwell  equation  (2-h)  through  the 
time  dependence  of  the  coupling  parameter  gi(t),  consistent  with  the  results  of  Ref.  7  and 
Ref.  13.  We  assume  this  dependence  to  be  exponential  in  time  with  a  decay  rate  or 
dephasing  rate  of  Fe  =  aF3. 

The  relation  of  the  dephasing  factor  to  the  individual,  slightly  different,  resonance 
energies  at  different  nuclei  (inhomogeneous  broadening)  is  discussed  in  Chapter  n  of  this 
report  In  that  chaper  we  also  show  how  relaxation  can  reduce  the  broadening  effect  and 
thereby  recover  the  SF  pulse.  In  die  present  paper  we  study  the  production  of  the  SF  pulse 
as  a  function  of  b  or  a,  described  by  Eqs.  (2-a)-(2-i). 
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The  results  are  shown  in  Figs.  10a  and  10b  where  the  natural  decay  was  assumed 
to  be  r  =  10^.  At  low  values  of  b  or  a  (when  Fe  «  Td”1)  there  is  no  noticeable 
effect  on  either  the  delay  time  or  the  intensity.  At  around  b,  a  =  10^  (F^  Fe  =  the 
pulse  delay  increases,  and  reaches  a  maximum,  while  the  intensity  starts  to  decrease.  -For 
b,  a  >  1(P  the  delay  time  decreases  but  the  intensity  drops  to  the  spontaneous  non- 
cooperative  value.  The  increase  in  Td  occurs  during  the  time  when  the  dephasing  effect  of 
F^  or  Fe  partially  neutralizes  the  correladons  built  up  by  the  electric  field  vacuum  state 
fluctuations,  and  consequently  more  time  is  required  to  build  up  the  SF  pulse.  The 
decrease  of  both  Xd  and  lo  indicates  the  complete  breakup  of  SF  past  a  certain  time,  and 
thus  only  the  initial  part  of  the  SF  peak  appears,  but  becomes  smaller  as  a  or  b  increases. 
The  only  difference  between  the  effects  of  a  and  b  on  the  SF  pulse  is  the  stronger  reduction 
in  the  peak  intensity  lo  and  the  more  pronounced  increase  in  the  delay  time  xd  due  to  the 
inhomogeneous  broadening  as  compared  to  the  homogeneous  broadening,  for  the  same 
amount  of  broadening  a  =  b. 

Some  insight  into  the  difference  between  homogeneous  and  inhomogeneous 
broadening  in  producing  an  effect  on  the  SF  emission  can  be  obtained  by  considering 
equations  (2-f)  through  (2-i).  Restricting  ourselves  to  one  transition  (3  to  1)  only,  after 
dropping  the  subscripts  we  have 

^  =  (2-i)' 

where  f(x,  t)  =  2N3(t)  +  N4(t)  -  Nq. 

Equation  (2-g)'  and  (2-i)'  can  be  combined  together,  assuming  an  exponential  form 

of  g(t), 

^=[g(.)f(x..)-i(r+r^+r,)R]E* 
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ISO 


Broadening  Parameter  a  or  b 


Figure  10(a).  SF  Pulse  Emission  Delay  Time  xp  as  a  Function  of  the 
Homogeneous  Broadening  Parameter  b  (solid  line)  and  the  Inhomogeneous 
Broadening  Parameter  a  (dashed  line).  Curves  (a),  (b)  and  (c)  were  obtained 
with  N  s  10^,  10*,  10*.  respectively,  (b)  SF  Pulse  Emission  Intensity  Iq 
as  a  Function  of  the  Broadening  Parameters  a  and  b. 
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We  see  from  the  above  that  the  inhomogeneous  broadening  not  only  contributes  an 
additional  term  to  the  homogeneous  effect  but  also  through  the  time  dependence  on  g(t) 
reduces  the  coupling.  This  enhanced  destructive  effect  is  reflected  in  our  Monte  Carlo 
results.  Hgure  8  exhibits  the  effect  of  the  two  broadening  mechanisms  tm  the  full  SF  pulse 
shape  when  a  and  b  take  on  the  values  200,  500,  and  1000.  In  Figs.  11(a),  (b),  and 
(c)  a  =  0  and  in  Figs,  (d),  (e)  and  (f)  b  =  0.  In  Fig.  1 1(a)  the  dashed  Une  represents,  for 
comparison,  the  pulse  with  no  broadening  a  =  b  =  0,  and  the  solid  line,  the  pulse  when  a  = 
0,  b  =  200.  As  expected,  from  Figs.  10(a)  and  (b)  the  pulse  intensity  decreases  more  with 
inhomogeneous  broadening  and  the  pulse  spreads  out  more  than  with  the  homogeneous 
broadoiing.  The  delay  tintt  is  also  affected  differently  at  these  values  of  the  parameters. 


5-2-91 -5 


Figure  il.  SF  Pulse  Emission  Intensity  lo  as  a  Function  of  the  Homogeneous 
Broadening  Parameter  a  and  the  Inhomogeneous  Broadening  Parameter  b. 
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E.  TEMPORAL  VARIATION  OF  FIELDS  INSIDE  THE  ACTIVE  MEDIUM 


The  solution  of  Eqs.  (2-a)-(2-h)  is  a  function  of  space  and  time.  We  have  so  far 
discussed  the  emission  of  the  SF  pulse  from  the  right  face  of  the  active  medium  as  a 
function  of  the  characteristics  of  the  medium.  This  is  what  would  be  measured  in  an 
experiment  It  is  interesting  to  make  use  of  the  infcKmation  about  the  activity  in  the  medium 
as  the  pulse  develops  and  propagates  through  the  medium.  We  are  considering  the 
geometry  shown  in  Fig.  12. 


Z 


Figure  12.  Schematic  of  the  Geometry  Used  in  the  Calculations  of  the  Internal 

Fields  in  the  Active  Region. 


The  fields  inside  the  medium  (o  <  x  <  1)  are  functions  of  x  and  t,  polarized  in  the  z 
direction.  The  fields  at  x  =  1  give  the  emitted  pulse  propagating  in  the  positive  x  direction. 
The  results  of  our  study  are  summarized  by  eight  runs  for  the  parameters  given  in  Table  2. 

The  result  of  the  eight  runs  are  shown  in  three  types  of  figures.  One  set  of  figures 
(type  1)  gives  the  temporal  variation  of  the  nonnalized  level  densities  N4,  with  initial  values 
N4(0)  =  1,  N3(t)  and  Ni(t)  with  initial  values  zero,  E+Ct)  and  I(t).  These  parameters  are 
plotted  in  Figs.  13-20  with  N4,  N3,  N],  E+  and  1  plotted  as  graphs  (a)-(f)  of  these  figures, 
respectively.  The  second  set  of  figures  (type  2)  gives  the  spatial  variation  of  E+(x),  R‘*‘(x) 
and  I(x)  at  selected  times.  The  electric  field  E***  for  the  different  runs  is  plotted  in 
Figs.  20(a)-28(a),  polarization  R+(x)  in  Figs.  20(b)-28(b)  and  the  intensity  in  Figs.  20(c)- 
28(c).  The  final  set  of  figures  (type  3)  gives  three  dimensional  plots  of  intensity  I(x,t),  as  a 
function  of  position  in  the  active  medium  and  time.  These  plots  are  given  in  Figs.  28-36 
for  runs  1  through  8,  respectively,  with  the  parameters  for  these  runs  given  in  Table  2. 


Table  2.  Parameter  Values  Characterizing  the  SF  Calculations  of  Runs  1*8 


Parameter 

Run  No. 

1 

2 

3 

4 

5 

6 

7 

8 

Y 

10 

10 

1 

BM 

10 

10 

1 

10 

r 

10-" 

10"^ 

1 

o 

10"“* 

10“^ 

1 

o 

10~^ 

10-3 

0 

0 

0 

0 

10 

5 

5 

5 

N 

,0= 

10® 

10® 

io" 

10® 

10® 

io" 

10" 

For  runs  with  p.  =  0  (1, 2, 3, 4)  the  spatial  variations  of  the  fields,  in  particular  the 
intensity,  is  monotonically  increasing  with  an  increasing  slope,  whereas  for  runs  with 
|i.  ^  0  (runs  S,  6)  the  electric  field  changes  sign;  consequently,  at  some  point  in  the  region 
the  intensity  goes  to  zero.  From  Eq.  (2-h),  which  gives  the  only  gradient  term,  me  can  see 
that  this  happens  when  the  build  up  in  the  polarization  and  the  electric  field  are  such  that 

X  X 

E(x)  =  g,jR(x)-|ljE(x).0  . 

O  O 

This  result  clearly  indicates  that  for  maximum  intensity  out  of  the  active  region  the 
length  of  the  region  has  to  be  carefully  selected  with  consideration  of  the  attenuation  p  in 
the  medium  and  the  coupling  constant  gi  between  the  nuclei  and  the  electromagnetic  Held. 
The  actual  position  of  the  null  in  the  intensity  is  not  an  obviously  simple  function  of  p.  and 
cannot  be  obtained  from  analyses,  such  as  the  triple  integral  approach,  which  average  over 
the  spatial  variable.  In  runs  7  and  8  there  is  a  minimum  but  not  a  complete  extinction  of  the 
field  at  some  x  value.  Further  analysis  of  the  transport  phenomena  is  required  to 
understand  the  full  significance  of  the  effect  of  propagation  in  the  active  region  and  the 
generation  of  the  SF  pulses  for  arbitrary  conditions. 


Figure  13.  SF  PuIm  Characteristics  for  Input  Parameters  of  Run  1  (see 
Table  2).  The  top  two  curves  show  the  variation  of  level  densities 
N4  and  Na,  calculated  at  x  s  /.  The  bottom  two  curves  give  the  level 
density  Ni  and  the  pulse  Intensity,  (y  s  10.  r  s  i(r^,  4  s  0.  N  s  i(^) 
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RUN  2 


Flgura  14.  SF  PulM  Characteristics  for  Input  Parameters  of  Run  2  (see 
Table  2).  The  top  two  curves  show  the  variation  of  level  densities 
Ns  and  Ns,  calculated  at  x  s  I.  The  bottom  two  curves  give  the  level 
densHy  Ni  and  the  pulse  Intensity,  (y  «  io,  r  s  p  s  0,  N  »  lO^) 
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FIgura  16.  SF  PulM  Characteristics  for  Input  Parameters  of  Run  4  (see 
Table  2).  The  top  two  curves  show  the  variation  of  level  densities 
N4  and  Ns,  calculated  at  x  s  /.  The  bottom  two  curves  give  the  level 
density  Ni  and  the  pulse  Intensity,  (y  r  0.1,  r  s  10-^,  p  s  0,  N  &  1(^) 
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Figur*  17.  SF  Pulse  Characteristics  for  input  Parameters  of  Run  5  (see 
Table  2).  The  top  two  curves  show  the  variation  of  level  densities 
Ns  and  Ns,  calculated  at  x  s  /.  The  bottom  two  curves  give  the  level 
density  Ni  and  the  pulse  Intensity,  (y  s  io,  r  s  itr^,  p  s  io,  N  s  lO^) 
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Figun  18.  8F  PulM  Charaet«rl«tlcs  for  Input  Paramotors  of  Run  6  (too 
T^Ha  2).  Tha  top  two  eurvaa  show  tha  variation  of  laval  dansitlaa 
N4  and  Ns,  caleulatad  at  x  ■  I.  Tha  bottom  two  eurvaa  glva  tha  laval 
danalty  Ni  and  tha  pulaa  Intanalty.  (7  «  10,  r  «  i(r^,  11  ■  s.  N  « 
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Figura  19.  SF  PulM  Characteristics  for  Input  Parameters  of  Run  7  (see 
Table  2).  The  top  two  curves  show  the  variation  of  level  densities 
Ns  and  Ns,  calculated  at  x  x  /.  The  bottom  two  curves  give  the  level 
density  Ni  and  the  pulse  imensKy.  (y  x  i,  r  x  io~2,  p  x  5,  N  x  lO^) 
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Figure  20.  SF  PulM  Charecttrlstics  for  Input  Paremotore  of  Run  8  (saa 
TatNa  2).  Tha  top  two  curvaa  aftow  ttia  variation  of  laval  danaltlaa 
N4  and  Na,  caleulatad  at  x  «  f.  Tha  bottom  two  curvaa  giva  tha  laval 
danaKy  Ni  and  tha  pulaa  intanaKy.  (y  «  10,  r  >  lO-^,  p  s  5,  N  «  10^) 


Figure  21.  Plot  of  the  Spatial  Variation  of  the  Polarization  R'*’(x),  Number 
Density  N3(x),  Field  Ni(x),  and  intensity  l(x)  at  Different  Times  Indicated  by 
Frame  Numbers.  The  times  corresponding  to  frame  numbers  are  given  by 
2.5  TsF  ^  (frame  number)  and  the  pulse  corresponding  to  these  data 
(Run  1)  as  shosm  In  Fig.  13.  {y  s  10.  r  s  1(r*,  p  s  0,  N  s  10^) 
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Figura  22.  Plot  of  '^ho  Spatial  Variation  of  tha  Polarization  R'^(x),  Number 
Denaity  N3(x),  Field  Ni(x),  and  Intenaity  l(x)  at  Differem  Timea  indicated  by 
Frame  Numbera.  The  timee  corresponding  to  frame  numbers  are  given  by 
3.75  tsF  X  (frame  number)  and  the  pulse  corresponding  to  these  data 
(Run  2)  as  shown  in  Fig.  14.  (y  s  10,  r  s  10~^,  p  s  0,  N  s  10*) 
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Figuni  23.  Plot  of  the  Spatial  Variation  of  the  Polarization  R'^(x),  Number 
Denalty  NaCx),  Field  Ni(x),  and  Intensity  l(x)  at  Different  Times  Indicated  by 
Frame  Numbers.  The  times  corresponding  to  frame  numbers  are  given  by 
3.75  tsF  *  (frame  number)  and  the  pulse  corresponding  to  these  data 
(Run  3)  as  shown  in  Fig.  15.  (y  s  i,  r  s  ^  s  o,  N  s  lo^) 


FIgura  24.  Plot  of  the  Spatial  Variation  of  the  Polarization  R'*’(x),  Number 
Density  N3(x),  Field  Ni(x},  and  Intensity  l(x)  at  Different  Times  Indicated  by 
Frame  Numbers.  The  times  corresponding  to  frame  numbers  are  given  by 
8.0  TsF  X  (frame  number)  and  the  pulse  corresponding  to  these  data 
(Run  4)  as  shown  In  Fig.  16.  (y  s  0.1,  r  s  10~^,  p  s  0,  N  s  10^) 


Figure  25.  Plot  of  the  Spatial  Variation  of  the  Polarization  R+(x).  Number 
Density  N3(x),  Field  Ni(x).  and  Intensity  l(x)  st  Differem  Times  Indicated  by 
Frame  Numbers.  The  times  corresponding  to  frsme  numb  'rs  are  given  by 
7.5  xsF  X  (frame  number)  and  the  pulse  corresponding  lo  these  data 
(Run  5)  as  shown  In  Fig.  17.  {y  *  10,  r  s  10-^,  p  s  10.  N  s  10®) 
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Figure  26.  Plot  of  th«  Spatial  Variation  of  tha  Polarization  R'^(x),  Number 
Denalty  N3(x),  Field  Ni(x),  and  Intensity  l(x)  at  Different  Times  Indicated  by 
Frame  Numbers.  The  times  corresponding  to  frame  numbers  are  given  by 
3.75  TsF  X  (frame  number)  and  the  pulse  corresponding  to  these  data 
(Run  6)  as  shown  In  Fig.  18.  (y  s  io,  r  s  u  s  5,  N  s  io?) 


Figure  27.  Plot  of  the  Spatial  Variation  of  the  Polarization  R-^(x).  Number 
Denalty  N3(x),  Field  Ni(x),  and  Intenelty  l(x)  at  Different  Times  Indicated  by 
Frame  Numbers.  The  times  corresponding  to  frame  numbers  are  given  by 
4.0  xsF  X  (frame  number)  and  the  pulse  corresponding  to  these  data 
(Run  7)  as  shown  in  Fig.  19.  (y  s  1,  r  s  lO"*,  p  s  5.  N  s  io3) 
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Figurt  28.  Plot  of  tho  Spatial  Variation  of  the  Polarization  R'*‘(x),  Number 
Denalty  N3(x),  Field  Nt(x),  and  intensity  l(x)  at  Different  Times  Indicated  by 
Frame  Numbers.  The  times  corresponding  to  frame  numbers  are  given  by 


5.0  tsF  *  (franid  number)  and  the  pulse  corresponding  to  these  data 
(Run  8)  as  shown  In  Fig.  20.  (y  s  10,  r  s  10~3,  u  s  5.  N  s  lO^) 


Figura  29.  intansHy  as  a  Function  of  Space  and  Time  for  Run  1. 

(Y  S  10,  r  *  10-^,  p  «  0,  N  B  103) 

The  time  axle  la  labeled  In  frame  numbers  (50  frames  s  50  tsp). 
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Figura  31.  Intonslty  as  a  Function  of  Spaca  and  Tima  for  Run  3. 
(y  .  1,  r  .  10-^  |i  «  0.  N  «  103) 

Tlia  tima  axis  Is  labalad  In  frama  numbars  (50  framas  «  75  tsf). 
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Figur*  32.  Intensity  m  a  Function  of  Spacs  and  Tima  for  Run  4. 

(y  c  0.1,  r  >  10-^,  4  s  0.  N  «  10*) 

Tlw  tima  axis  la  labalad  in  frama  numbars  (50  framat  •  100  tsp). 
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Figure  33.  Intensity  as  a  Function  of  Space  and  Time  for  Run  5. 

(T  *  10,  r  S  10-^,  |i  B  10.  N  S  10*) 

The  time  axis  Is  labeled  In  frame  numbers  (50  frames  s  150  tsp). 
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Figurt  35.  Intansity  as  a  Function  of  Spaca  and  Tima  for  Run  7. 
(Y  r  1,  r  s  10-2,  |i  s  5,  N  s  103) 

Tha  tima  axis  is  labalad  In  frama  numbars  (50  framas  a  200  xsf). 
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F.  CONCLUSIONS 


This  paper  considers  a  generalized  version  of  the  Haake-Reibold  model  of  SF  in  a 
multilevel  system  of  emitters,  having  up  to  five  energy  levels.  The  generalized  version 
takes  into  account  finite  pumping  rates,  spatial  attenuation  of  the  electromagnetic  field, 
dephasing,  and  population  depletion  due  to  processes  such  as  internal  conversion  and 
spontaneous,  but  noncollective,  transverse  emission. 

In  this  paper  we  studied  the  dependence  of  emitted  pulse  sh^s  on  the  statistics  and 
time  dependence  of  the  noise  sources.  Monto  Carlo  results  obtained  by  averaging  25  to 
100  pulses,  which  seem  to  be  adequate  samples,  also  depend  upon  the  noise  source 
variance. 

The  spatial  attenuation  parameter  p  reduces  the  pulse  peak  and  shifts  it  to  longer 
times.  The  peak  will  decrease  gradually  but  will  not  disappear  if  no  dephasing  mechanisms 
that  reduce  the  dipole  built  up  by  the  noise  source  exist  and  no  population  depletion 
mechanisms  destroy  the  inversion. 

Small  values  of  the  dephasing  factor  have  no  effect  on  the  emitted  pulse  delay  time 
or  intensity,  but  when  the  dephasing  rate  is  comparable  to  the  SF  decay  1/tSR  rate  the 
delay  time  starts  to  increase  with  increasing  values  of  the  dephasing  factor  until  it  reaches  a 
maximum,  after  which  it  decreases  rapidly. 

Because  of  spatial  attenuation,  for  maximum  intensity  of  the  pulse  radiated  out  of 
the  active  region  the  choice  of  the  region’s  length  should  be  governed  by  the  values  of  the 
attenuation  factor  p  and  the  coupling  constant  gi  associated  with  coupling  between  the 
emitters  and  the  electromagnetic  field.  The  optimum  length  will  depend  on  details  of  the 
propagating  electromagnetic  field  inside  the  active  region. 
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n.  AN  ANALYSIS  OF  THE  EFFECTS  OF  INHOMOGENEOUS 
BROADENING  AND  TIME-DEPENDENT  INTERACTIONS 
ON  SUPERFLUORESCENT  EMISSION 


A.  INTRODUCTION 

If  a  "gamma-ray  laser"  is  developed  it  may  emit  in  a  superradiant  (SR)  or 
superfluorescent  (SF)  iiKxle  instead  of  a  stimulated  emission  (SE)  mode  (Refs.  4,  14). 
There  are  some  characteristic  features  of  superfluorescent  decay  which  may  provide  ways 
of  overcoming  specific  nuclear  problems  (Ref.  3)  which  at  present  prevent  CW  or  pulsed 
laser  (SE)  action.  In  another  paper  (Ref.  3)  we  have  investigated  the  possibility  of  nuclear 
superfluQiesoence  and  described  the  intenelation^p  of  the  pertinent  nuclear  and  solid-state 
parameters  that  govern  the  experimental  realization  of  this  phenonoenon.  We  also 
contrasted  the  nuclear  SF  problem  with  the  atcmic  and  molecular  SF  problem  and  described 
several  special  issues  that  differentiate  the  two. 

The  differences  are  due  to  (a)  the  shorter  wavelength  of  the  emitted  nuclear 
radiation,  (b)  the  subsequent  large  recoil  energy  compared  to  the  natural  linewidth,  (c)  the 
generally  higher  multipolarity  of  the  nuclear  radiation,  (d)  the  common  occurrence  of 
completing  transitions  in  nuclear  processes,  (e)  the  relatively  high  electronic  attenuation  for 
this  nuclear  radiation  con^ared  to  the  resonant  cross-section,  (f)  the  difficulty  of  pumping, 
requiring  finite  pumping  times  for  nuclear  systems,  (g)  the  difficulty  of  providing  desired 
geometrical  shapes  for  the  active  medium,  and  (h)  inhomogeneous  broadening  which  can 
destix^  the  resonant  condition  when  narrow  nuclear  lines  are  involved. 

All  of  these  issues  are  of  secondary  interest  in  atomic  SF  or  SR  studies  but  they 
present  critical  problems  to  the  development  of  nuclear  SF  or  SR,  as  well  as  lasing,  or 
amplified  spontaneous  emission  (ASE). 

In  a  recent  paper  we  discussed  SF  at  wavelengths  shorter  than  the  interatomic 
spacing  (Ref.  IS),  the  effect  of  competing  transitions  and  electronic  attenuation  on  the  SF 
pulse  generation,  and  analyzed  the  effect  of  firtite  pumping  times  on  the  generation  of  the 
SF  pulse  (Ref.  3).  In  this  paper  we  discuss  how  inhomogeneous  broadening  affects  the 


generation  of  a  SF  pulse  and  how  the  effect  can  be  reduced  through  electronic  relaxation  or 
time-dependent  hypeifine  interactions. 

In  this  paper  we  discuss  the  calculation  of  the  emission  line  shape  under  dynamic 
conditions  (the  existence  of  time-dependent  hyperfine  interactions)  and  then  generalize  the 
procedure  to  a  system  of  many  nuclei.  We  calculate  the  emission  pulse  shape  for  a  system 
of  nuclei  in  slightly  different  environments,  which  causes  inhomogeneous  broadening,  and 
then  study  the  effect  of  time-dependent  perturbations  on  the  line  shape  and  pulse 
generation. 

B.  SINGLE  NUCLEUS  EMISSION  LINESHAPE 

Consider  a  two-state  nuclear  system  with  energy  levels  a  and  X  as  shown  below 
connected  by  an  electromagnetic  transition: 

- i - X  . 


The  probability  for  the  emission  of  a  single  photon  with  energy  h  (O,  wave  vector 
^  and  polarization  p  is  given  by  the  Wigner-Weisskopf  formula 


Pxa  = 


l<al  V^^^X>P 


(6) 


where  F  is  the  natural  line  width  for  the  transition  from  state  a  with  energy  Ea  to  state  X 

(+)  — ^ 

with  energy  E^,  and  V  (kcap)  is  the  Hamiltonian  for  the  emission  of  a  photon  with 
parameters  It,  ca,  p. 

At  this  point  we  can  make  some  modification  to  Eq.  (6)  which  will  allow  us  to  use 
well-known  results  and  convenient  formalisms  to  deal  with  multiatom  (or  nuclear) 
emissions  in  a  time-dependent  environment  The  first  step  is  to  rewrite  Eq.  (6)  in  terms  of 
amplitudes  and  the  second  is  to  transform  it  into  the  time  domain.  Equation  (6)  can  be 
rewritten  in  terms  of  the  amplitude  Axa  for  the  process  of  interest  as  follows: 


where 


^Xa~ 


<Xlv^'Na> 


<b  +  E„-E,1- 


Now  using 


f  -i(0»+Eo— Ej.)t'— ■r-t' 


y -  =  i  fdt' 

[w  +  Ea-Ej^j-i  j  o 
we  take  the  second  step  and  get 

T  .  (+).  -i(«>+Ea-E^)r-f  f 

A^=  i j  dt'<XiV  |a>e 

O 

=  jdt’<3lle"*  V^‘'^e-‘"*la>e““  2* 

O 

=  J  dt'<X\y^*\o  Ia>e 


-ieDfX  f 


(8) 


(9) 


(10) 


where 


since 


f/IX>  =  Ej^|X>and//la>  =  Eola>  . 


(11) 


To  obtain  a  lineshape  for  a  single  resonance  (transition  a 
wnitring  atom  at  position  i  we  write 

P.^(op)  =  A*  .A, 
i '  Xoi  Xa 

~  «  i<o(i'-t>f(f-t-) 

=  J  df  f  dt"e  ^ 

o  o 


X)  from  a  single 
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<X  |V^'\t’)vV)la>  . 


(12) 


Generalizing  to  multiple  level  systems  with  allowance  for  transitions  between  levels 
we  write 


oe  90 

P.  (k^mp)  =  J  dt’  J  dt"  e 


(f+t")  (.)  (+) 

2  <VV)V  (t’)>  .  (13) 


where 


V%)  =  e“v®e“'' 


and  the  bracket  pair  notation  means 

< A >  =  tr <pA >  =  ^  <XlplX’><X’lA|X>  , 
X.  X' 


so  that 

(-)  (+)  (-)  (+) 

<V  (t")V  (t')>  =  <V  V  (t*-t")>  . 

Introducing  a  change  of  variables  t  =  t'  —  t"  and  integrating  with  respect  to  t”  one 
gets  a  convenient  form  for  the  emission  lineshape  assuming  V(t)  is  a  stationary  process 

90  .  r 

^  0  f  (-)  (+) 

P.  (^(op)  Re  I  dx  e  ^  <  V  V  (x)  >  .  (14) 

O 

This  form  is  a  starting  point  for  a  description  of  dynamic  effects  (relaxation)  in 
Mossbauer  emission  and  absorption  experiironts.  Such  experiments  are  usually  modeled 
assuming  nuclei  emitting  noncooperatively.  An  example  of  line  shape  calculations  for  a 
specific  relaxation  mechanism  will  be  presented  in  Section  F  below. 


C.  EMISSION  FROM  A  SYSTEM  OF  N  NUCLEI 


For  N  atoms  using  Eqs.  (1 1)  and  (10)  we  get 

12 


Pxa(“)=  aL  = 


‘Xoi 


f  I  ifi.l'  {-)  -iHX , 

2  J  dt’<Xj  le  *  V  e  ‘  ' 


a.  >e 


r 

t* 

2 


1  o 


(15) 


Subscript  i  is  used  to  indicate  that  even  though  the  states  are  the  same,  the 
energies  can  be  slightly  different  for  different  nuclei,  resulting  in  inhomogeneous 


faroadening  of  the  line  shape  for  the  total  system  (Ref.  16).  For  a  stationary  case  we  use 
//.  |aj>  =  laj>toget 


^  J  df<A.j 
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=  J  dt’  J  dt"  <  \  I I  a  ><  aj  i  Xj  > 

i  o  o 

^  - 1")  +  ko(f  - 1")  -y  (t‘  + 1") 

+  ^  J  df  J  dt"  <X.J  V‘^laxa.1 

i.  i  o  o 

r 

MM*  -  ieOjt"  +  in  (t’  - 1")  (t'  + 1") 

X  e  ‘  ^ 

ii  A  A  A 

=  II«0.©.)+  I  I(0);(D..tD.)  =  I  ,+ I  ,  (16) 

i  ‘  ij^ti  *  J  N  N 

The  first  term  in  Eq.  (16)  gives  the  lineshape  for  spontaneously  emitted  photons 

fiom  the  N  nuclei  acting  independently  and  is  equal  to  ^  P.  ^cop)  of  Eq.  (14).  The 

i 

second  term  gives  the  lineshape  for  the  collective  emission.  Each  term  in  the  sum  is  a 
Lorentzian  with  slightly  different  resonance  energy  (Di.  The  sum  gives  an 
inhomogeneously  broadened  line  which  usually  is  assumed  to  have  a  Gaussian  or 
Lorentzian  envelope  and  with  slight  modifications  is  appropriate  to  describe  both 
Mdssbauer  absorption  and  emission  experiments  (Refs.  17,  18).  In  the  special  case  of  a 
single  nucleus,  Eq.  1  is  the  starting  point  for  Blume's  discussion  of  the  effect  of  time 
dependent  hyperfine  interactions  on  the  Mdssbauer  lineshape  (Ref.  19).  In  the  case  of  N 
nuclei  with  appropriate  normalization,  it  represents  the  source  lineshape  and  also  the 
absorption  cross-section  in  the  calculation  of  the  effect  of  inhomogeneous  broadening  on 
Mdssbauer  experiments. 
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D.  CALCULATION  OF  THE  EMISSION  PULSE  SHAPE 


For  discussing  SF  emission  we  are  primarily  interested  in  the  total  intensity  as  a 
function  of  time  rather  than  the  lineshape  or  photon  energy  distribution.  After  inversion  is 
accomplished,  the  delay  time  xd  tuid  the  emission  time  Tsf>  both  shorter  than  the  natural 
lifetime  of  a  single  atom  or  nucleus,  characterize  the  SF  process. 

With  this  in  mind  we  consider  the  second  term  of  Eq.  (16) 

I  2(o>)=  Z  1(0);  (D.,  0).)  =  S  J  dt*  J  dt"<X.  Iv^^la.  > 

PT  '  ■*  o  o  *  * 

r 

i  co.t*  - 1  m.t"  + 1(0  (t’-t")  -  y  (f  + 1") 

x<a.V  ‘  ^  (17) 

1  J 


We  integrate  I  ,  (o))  over  O),  which  produces  a  delta  function  in  (t'  - 1")  and  allows 

tr 

us  easily  to  compute  the  integral  over  t",  for  the  total  emission  intensity. 
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.(-) 


i  .^1  i  df<X  \w  la.xalv 
PT  i  o  *  •  J 


,(+) 


X.  >  e 
J 


i((0.-  (0^)  t'  -  r  t’ 


Accordingly,  for  the  emission  intensity  as  a  fiiiKnion  of  time  we  assume 

I  ,(t)  =  y'  <X.|V‘^|a.><a.|V''^|X.> 

PT  *  1  j  j  e 

i,j^i 


(18) 


As  can  be  seen  from  Eq.  (18)  the  effect  of  inhomogeneous  broadening  from  the 
difference  o)i  -  ooj  in  the  exponential  is  to  introduce  a  phase  change  between  different 
radiators.  This  lowers  the  intensity  of  the  emitted  SF  pulse.  In  Mbssbauer  spectroscopy 
the  effect  of  inhomogeneous  broadening  is  thought  of  as  a  reduction  in  overlap  of  source 
and  absorber  lines  due  to  shifting  in  the  individual  resonances  by  detuning  or  different 
resonance  energies  at  different  sites.  This  is  discussed  in  Ref.  17. 

In  discussions  of  lasing  the  effect  of  inhomogeneous  broadening  is  used  to  describe 
the  reduction  of  the  stimulation  cross-section  and  consequently  the  loss  of  gain  in  the 
operation  of  a  laser.  However,  in  the  SF  phenomena,  as  pointed  out  by  Eberiy  (Ref.  13)  it 
is  the  phase  change  introduced  by  the  inhomogeneous  broadening  that  controls  the 
emission  and  not  the  reduction  of  overlap  of  resonance  lines.  For  treating  SF  emission  it  is 


the  phase  change  introduced  by  different  resonant  energies  0)i  at  different  sites  i  that 
describes  the  effect  of  inhomogeneous  broadening  on  the  emission. 

At  this  point  it  is  instructive  to  compare  Eberly's  derivation  of  the  SF  emission  as 
effected  by  the  inhonoogeneous  broadening  and  our  results  given  in  Eq.  (16).  Eberly 
shows  that  the  intensity  of  emitted  radiation  as  a  fiuiction  of  time  t  and  the  wave  vector  k  is 
given  by 

i(k;t)=i„fl?)  <x  iiCr^ix  > 

ni=  1 

+  y  X  <x 

e  m^e 


c 


(19) 


where  in  Eberly's  notations  R^,  are  the  usual  raising  and  lowering  operators  for  the 
mdi  nucleus  or  atom,  I  x  >  is  the  atomic  (or  nuclear)  state  in  the  interaction  picture  so  that 

I X  >  =  U(T,  0)  I  Vh  >  with  U(J,  0)  as  the  unitary  operator  and  I  Vh  >  the  Heisenberg 
state.  The  detuning  frequencies  Ym  arc  obtained  from  Ym  *  Vm  -  co  where  Vm  is  the 
transition  frequency. 

The  first  sum  in  Eq.  (19)  gives  the  sin^e  nucleus  emission  and  is  related  to  the  first 
part  of  Eq.  (16).  The  second  sum  in  Eq.  (19)  gives  the  SF  emission  with  the  characteristic 
focused  angular  distribution  and  shortened  emission  time  following  an  initial  delay.  We 
will  now  concentrate  on  the  second  sum  and  consider  the  effect  of  inhomogeneous 
broadening. 


Eberly  shows  that 


I  I 

l  m*l 


)» 
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-  N 


n-Y 


and  the  sum  squared  term  can  be  written  as 


|2 


I 


i(D  t 


=  N^ 


r  ica  t1 
|c  "  I 

‘  ^ave 


(20) 


where  is  the  inhomogeneously  broadened  lineshape  function  normalized  to  unity.  In 
the  case  when  g(Y)  can  be  modeled  by  a  Lorentzian  lineshape,  the  Fourier  transform  in 
Eq.  (20)  leads  to  exponential  damping,  and  the  SF  intensity  term,  the  part  of  Eq.  (16), 
becomes 

I  ,(t)  =  brl<V>  1  e  *  (21) 


*1^  is  the  inverse  of  the  inhomogeneously  broadened  line  width.  For  fast  pumping, 

essentially  completed  near  t  =  0  (such  that  tpump  « ‘tSF  « *^2^  inhomogeneous 
broadening  has  little  effect  on  the  emission  bectuise  the  nuclei  have  not  had  time  to  dephase. 

Only  after  a  relaxation  time  does  the  phase  difference  between  different  nuclei  i,  j, 
(Oii  -  09j)t,  become  large  enough  to  affect  the  emission  or  to  prevent  the  SF  pulse  from 
occurring. 

Thus  inhomogeneous  broadening  affects  SF  pulse  emission  differently  from  the 
way  it  affeas  lasing  or  amplified  spontaneous  emission  (ASE).  The  difference  between 
lasing  and  superfluorescence  can  be  understood  in  terms  of  the  different  times  during  the 
evolution  of  a  system  that  they  dominate  the  emission  and  the  kinds  of  processes  that  they 
resemble.  SF  is  like  an  emission  process  that  is  cooperative  in  nature,  whose  characteristic 
is  described  by  a  delay  time  to  and  the  SF  emission  time  tsp-  Lasing  is  a  process  like 
absorption  that  is  characterized  by  a  cross-section,  the  stimulation  cross-section  0$.  SF 
(under  appropriate  conditions)  can  occur  early  in  the  evolution  of  a  system  before  complete 
dqjhasing,  while  lasing  can  occur  even  after  dephasing  if  the  gain  coefficient  is  greata  than 
zero  (i.e.,  amplification  is  possible).  This  difference  in  the  effect  of  inhomogeneous 
broadening  cm  lasing  (ASE)  and  SF  has  implications  for  y-ray  laser  feasibility. 

Both  of  the  processes  have  to  be  considered  in  a  complete  study  of  f-xvj  lasers. 
Candidate  selection  should  be  made  using  the  criteria  for  both  ASE  and  SF  pulse 
generation.  If  inversion  can  be  accomplished  quickly  and  the  inversion  density  is  high 
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# 


*  J 

enough  to  initiate  SF  befcwe  dephasing  sets  in,  xsp  «  ,  SF  will  occur.  Otherwise, 

if  the  gain  is  high  enough,  ASE  will  occur.  The  transition  region  has  been  observed 
expoimentally  (Ref.  20)  in  atomic  systems. 

For  lasing  to  occur  it  is  necessary  that  the  photon  gain  per  unit  length,  K,  be 
greater  than  zero  (Refs.  21, 22). 

K  =  <yRAn*-4  . 
and 

2it(l  +  o)(l  +  a)  ■ 

with  X  the  wavelength  of  the  photon  emitted  in  the  resonant  transition,  f  the  recoilless 
fractitMi,  gN  the  coupling  cmistant,  a  the  internal  conversion  coefficient,  ^  the  branching 
ratio,  and  p  is  the  linear  attenuation  coefficient  Inhomogeneous  broadening  is  included 
throu^  die  term  (1  +  a)  and  the  parameter  a  can  be  10^  or  higher  for  isomeric  nuclei  of 
interest 

Thus  when  considering  lasing  or  ASE  the  stimulation  cross-secdon  is  assumed  to 
be  reduced  by  a  dme-indqiendent  factor  .  —  ,  while  in  SF  the  reduction  of  the  intensity  is 

1  T"  ft 

characterized  by  a  time-dependent  phase  factor  which  reflects  the  dynamic  nature  of  the 
process  and  the  evoluticm  of  the  inverted  state.  The  dephasing  process  can  be  described  in 
terms  of  the  dephasing  time  T2*  as  in  Eq.  (21). 

E.  TIME-DEPENDENT  PERTURBATIONS 

If  the  system  under  consideration  is  subjected  to  time-dependent  perturbing  effects, 
Eq.  (14)  can  still  be  used  to  calculate  the  emission  and  absorption  lineshapes  for  resonators 
acting  independently.  Assuming  a  time-dqjendent  field  acting  on  the  nuclear  spins  we  can 
write  H(t)  =  //q  +  HIj  f(t),where  f(t)  is  a  random  process.  The  lineshape  is  given  by 

r 

P(ka)p)=:|rRe/  dte  ^  . 


(22) 


(23) 


where  the  bar  indicates  an  average  over  all  functions  f(t). 


Let's  assume  that  f(t)  =  ±1,  selected  at  random  for  different  t.  and  H  =  Hq  =  Mogo^ 
for  the  ground  state  and  H  =  Hi  =  4igi/t  for  the  excited  state,  with  p-j,  gj  the  respective 
magnetic  moments  and  geomagnetic  ratio,  and  h  is  the  strength  of  the  field,  and  let  IT  =1 
for  convenience.  Blume  (Ref.  19)  has  computed  the  indicated  average  to  be 


(t)>  =  e 


I  I  <Lm„|  m  >1^ 


The  average  on  the  right-hand  side  of  (24)  can  be  calculated  as  follows: 


cos  X  Q  t  + 


■^sinQtjc-Oi 


(25) 


where 

(goino-g,m,]^ 

Q* 

In  the  above,  Q  is  the  relaxation  rate  or  probability  of  jumping  or  flipping  between 
values  of  f(t)  =  +1  and  -1,  with  equal  probabilities  of  finding  the  system  in  either  +1 
or  -1.  For  a  more  general  case  assuming  different  probabilities  of  finding  the  system  in 
state  f(t)  s  +1  and  —1,  Balko  and  Hoy  (Ref.  23)  have  calculated  the  Mossbauer  lineshape  to 
be 

(oA^  -  Air©  +  aA-)]-(i  +  a)  i  Py] 

!(©)=  - - - L__ - J  ^  (27) 

A  V  +  3  ryj  +  [r©  +  «(A-"  “  oA')] 

where 

A'^=©  +  5,  A"*-©  +  5,  y  =  (a+ l)Q+f^rl 


a  =  (T/Q)  +  k  . 


To  detennine  the  effect  of  time-dependent  perturbations  on  the  SF  pulse  we 
ggneraliTft  Eq.  (18)  to  multiple  levels  by  introducing  the  density  matrix  formalism  as  in  the 
discussion  following  Eq.  (13). 


Then 
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The  matrix  elements  can  be  rewritten  using  Eq.  (11)  to  show  explicitly  the  time 
dependence,  and  the  stochastic  average  can  be  taken  as  in  Eq.  (23),  using  the  Hamiltonian 
ofEq.  (22). 


The  matrix  elements  are 
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We  have  used  the  condition  that  V(^)  connects  only  the  ground  and  excited  states  and  our 
time-dependent  Hamiltonian  //(t)  only  changes  the  en^gy  of  a  state  but  doesn’t  induce 
transitions  between  states. 


Integrating  over  od,  since  we  are  interested  in  all  energies,  we  get  a  5(t'  - 1")  from 
the  e^  -  n  Performing  the  integral  over  t",  taking  the  average  over  initial  states  and  the 
sum  over  final  states,  and  concentrating  our  interest  on  a  single  nuclear  transition,  thus 
dropping  die  sum  over  a  and  X,  we  get 

I  ,(0=  Z  <X.  Iv.  (0)  la. ><a.  Iv.  IX. >e'^ 

i.>i  *  ‘  ‘  J  J  J 


where  on  the  right-hand  side  the  bar  over  the  two  factors  in  parentheses  indicates  the 
average  value. 

The  random  variation  of  f(t')  at  nucleus  i  is  independent  of  the  variation  of  f(0  at 
nucleus  j.  Generally,  one  has  to  take  into  account  at  the  same  time, 

(a)  the  independent  fluctuations  at  nuclei  i  and  j,  and 

(b)  the  variation  in  the  two  sites  of  the  parameter  a  over  an  appropriate 
distribution. 


The  stochastic  averages  can  then  be  taken  as  in  Eq.  (24)  using  the  Hamiltonian  of 
Eq.  (22).  We  get 
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G(Q,  a.  t)  =  [cos  X  Gt  sin  X  Gtje-' 


■fe-I 


a  =  |gl^^l“go^^| 

as  expected  for  Q  »  a,  G(t)  •>  1,  and  the  inhomogeneous  broadening  and  phase 

destruction  of  the  superradiant  pulse  disappears,  while  for  G  «  a,  G(t) »  y  (e“**  +  e'““) 

and  die  maTimuin  phase  destruction  occurs  if  G(G,  (X,  t)  is  calculated  for  a  distribution  of 
a's  given  by  g(a)  which  represents  the  appropriate  inhomogeneous  broadening. 

If  we  assume  die  distributicm  is  gaussian  widi  variance  then 


and  die  average  value 


?j(G,  o,  t)  =  J  7(0,  a)  G(G,  a,  t)  do 
0 


should  replace  G(G,  a,  t)  in  the  expression  Eq.  (30)  to  obtain  11^2(0. 

In  the  following  three  figures  we  show  some  examples  of  the  effectiveness  of 
inhcnnogeneous  broadening  dephasing  reduction  factor  G(G,  a,  x)  as  it  depends  on  G 
and  a.  In  all  cases  F  =  1  s~i  is  assumed. 

In  Fig.  37  we  show  the  function  G(G,  g,  t)  for  parameters  G  =  0,  a  »  1, 10,  and 
1(X)  as  a  function  of  time.  The  horizontal  line  labeled  a -0  shows  no  reduction  in  the  SF 
pulse  intensity  because  there  is  no  inhomogeneous  broadening.  The  other  three  curves 
show  the  decrease  in  SF  pulse  intensity  when  different  dephasing  rates  due  to  different 
amounts  of  inhomogeneous  broadening  are  present  A  SF  pulse  with  a  delay  time  of 
1.6  sec,  which  would  appear  at  the  positions  of  the  arrow  labled  0  when  a  =  0,  would  be 
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reduced  to  about  1/2  of  this  value  when  a  -  1.  A  pulse  appearing  earlier  (at  arrow  2) 
would  suffer  less  of  a  reduction.  In  Fig.  38  we  show  a  plot  of  these  results  with  the  time 
axis  on  a  logarithmic  scale.  For  the  three  g  values  the  drop-off  is  clearly  indicated  at 
different  times.  Figures  39  and  40  show  the  reduction  factor  for  g  =  1(P  and  10^, 
respectively,  and  different  values  of  the  relaxation  rate  il.  The  curves  decay  rapidly  for 
small  values  of  and  slow  down  for  higher  values  approaching  the  limit  of  1.0  as  Q 
approaches  infinity.  The  results  of  Figs.  39  and  40  indicate  that  time-dependent  hyperfine 
interactions  can  reduce  the  effect  of  inhomogeneous  broadening.  As  the  relaxation  rate  £2 
increases,  the  SF  intensity  also  increases,  and  in  the  limit  of  £2  »  g  the  SF  pulse  intensity 
approaches  the  value  one  obtains  when  g=0.  This  SF  pulse  intensity  is  determined  by  the 
number  of  cooperating  nuclei  and  the  intensity  of  nuclear  properties.  Thus,  if  a  candidate 
exists  that  when  inverted  would  produce  SF  were  it  not  for  inhomogeneous  broadening, 
our  results  show  that  it  is  still  possible  to  get  a  SF  emission  if  a  fast  enough  relaxation  of 
the  levels  producing  the  inhomogeneous  txoadening  can  be  realized. 


Figure  37.  Reduction  of  SF  Pulse  Intensity  by  Inhomogeneous 
Broadening  with  Variance  g^ 
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Figure  38.  SemFlog  Plot  of  the  Results  Shown  In  Figure  37. 
The  dropoff  times  for  different  values  of  a  are  clearly  indicated. 
For  (a)  o  8  1,  (b)  o  s  io  and  (c)  a  s  loo. 


Figure  39.  Recovery  of  SF  Pulse  intensity  by  Relaxation  In  the  Presence  of 
Inhomogeneous  Broadening  with  o  sio^.  The  relaxation  rates  are 
(a)  Q  8  10,  (b)  £2  8  10^,  (c)  and  (d)  Q  8  10^ 
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Figure  40.  Recovery  of  SF  Pulse  Intensity  by  Relaxation  In  the  Presence  of 
Inhomogeneous  Broadening  with  o  s  10*.  The  relaxation  rates  are 
(a)  Q  s  10.  (b)  Q  s  10^,  (c)  n  :::  10*.  (d)  £2  s  10^  and  (e)  Q  =  10^. 

F.  CONCLUSIONS 

Our  calculations  show  that  the  effect  of  inhtnnogeneous  broadening  on  SF  enussitxi 
is  qualitatively  different  from  that  on  ASE  or  lasing.  The  effect  on  SF  is  to  cause 
dephasing  of  the  different  nuclear  dipoles  that  are  forming  correlations.  This  dephasing  is 
time  dependent  and  has  a  time  scale  that  goes  inversely  with  the  inhomogeneously 
broadened  line.  For  times  less  than  this,  SF  emission  is  less  affected  by  this  broadening 
than  ASE  or  lasing.  For  the  latter,  the  reduction  is  time  independent  and  the  stimulation 

cross-section  is  reduced  by  a  factor  of  —  where  a  =  2V2//t2  ^  a  is  the  variance  of  the 

1  +a  1 

distributitm  characterizing  the  inhomogeneous  broadening  and  F  is  the  natural  line  width. 

We  also  show  that  time-dependent  hyperfrne  interactions  of  the  nucleus  with  the 
electrons  can  induce  a  reduction  of  the  dephasing  and  an  increase  in  the  probability  of  SF 
over  a  large  time  range.  These  electronic  relaxation  effects  are  well  known  for  their 
lineshape  modification  as  observed  in  Mdssbauer  and  nuclear  magnetic  resonance  (NMR) 
experiments.  For  SF  they  provide  a  means  of  reducing  the  dephasing  due  to 
inhomogeneous  broadening  and  thus  may  be  useful  in  overcoming  one  of  the  more  difficult 
obstacles  to  nuclear  SF. 


at 
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m.  COHERENT  EFFECTS  AND  INHOMOGENEOUS 

BROADENING 


A.  INTRODUCTION 

In  a  recent  publication  (Ref.  24)  we  derived  an  expression  for  the  inhomogeneously 
broadened  lineshape  expected  in  a  Mdssbauer  absorption  experiment  Assuming  that  the 
nuclei  interacted  incoherently,  the  expression  for  absorption  was  given  as  a  product  of 
intensity  of  incoming  beam  and  the  single  nucleus  absorption  cross-section.  We 
summarize  the  results  here  and  then  derive  an  expression  assuming  coherent  processes  for 
scattering  of  emissions  from  a  multinuclear  source  by  a  multinuclear  absorber.  The 
calculation  presented  here  was  done  to  check  the  claims  presented  at  a  recent  ultrashort 
wavelength  laser  review  meeting  that  homogeneous  broadening  equally  and  intentionally 
produced  in  both  source  and  absorber  lines  could  produce  an  enhanced  Mdssbauer  effect 

The  experimental  condition  assumed  for  both  calculations  is  shown  in  Fig.  41(a), 
and  the  lineshapes  for  both  the  inhomogeneously  broadened  source  and  absorber  are 
depicted  in  Fig.  41(b). 

B.  INCOHERENT  APPROACH 

The  detected  loss  in  the  total  photon  count.  Id,  because  of  the  presence  of  the 
absorber  is  given  by 

*2  +«> 

Ij,=  J  dt  j I (E, r, S) n a(E, n dE  (34) 

tj  ^ 

where  n  is  the  number  of  resonant  nuclei  in  the  beam  per  cm^,  and  the  measurement  is 
perforated  during  the  time  interval  t2  —  tj.  I  (E,  F,  S)  is  the  source  radiation  lineshape  and 
g(E,  F)  is  the  resonant  cross-section  of  the  nuclei  in  the  absorber. 
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Figure  41.  (a)  A  schematic  showing  the  MOssbauer  transmission  experimental 
geometry,  (b)  A  graphic  representation  of  the  source  and  absorber  llneshapes. 
The  vertical  lines  represent  some  of  the  natural  lines  of  the  different  nuclei 
In  the  source  and  absorber  whose  envelopes  are  the  Inhomogeneously 
broadened  lines.  S  represents  the  doppler  shift  of  the  source 
and  ei  and  ej  the  resonance  energies  of  nuclei  I  and  J 
In  the  source  and  absorber,  respectively. 
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The  source  lineshq)e  is  given  by 


I(E,r.S) 


(e,) 


frl 

lo 

f. 

[E-(Eo-ej)-sf  +  (r/2)^ 


(35) 


where  Iq  is  the  flux  of  photons  from  the  source  in  cm“^  s  ^  in  the  direction  of  the  detector, 
Eo  is  the  resonant  energy  of  an  unperturbed  nucleus,  Ej  is  the  perturbation  of  the  energy  of 
nucleus  i  due  to  the  environment,  S  is  the  Doppler  shift  due  to  the  relative  motion  of  the 
source  and  absmber  nuclei  and  is  given  by  S  =  v/c  E,  with  velocity  v  where  c  is  the  speed 
of  light,  r  is  the  natural  linewidth^  and  f^  is  the  recoil-free  fnu:tion  in  the  source.  It  is 
assumed  in  Eq.  (34)  that  the  integration  time  is  short  enough  that  Iq  can  be  taken  as 
constant 


The  distribution  of  nuclear  environments  is  assumed  to  be  normal  with  width  A  so 
that  the  probability  of  finding  a  nucleus  with  energy  Ei  displaced  from  the  resonance  energy 
Eq  in  a  width  AEj  is  given  by 


P,(e,)--t 


Jn  A 


Ae, 


(36) 


Similarly,  for  the  effective  cross-section  of  the  absorber  we  have 

-iw 


o(EJl 


(E-(E„-£i)f  +  ^)’ 


Where  Gq  is  the  cross-section  on  resonance  and  is  given  by 

2  21+1 


1 


_  X, 

^0  “  2jt  2r+r  TTct 

8 


(37) 


4  and  Ig  are  the  nuclear  spin  quantum  numbers  of  the  excited  and  ground  states, 
respectively,  a  is  the  internal  conversion  coefficient  of  the  transition,  X  is  the  nominal 
wavelength  of  the  gamma-ray,  Fy  is  the  photon  absorption  partial  width. 


^  r  is  die  total  level  width  which  is  a  sum  of  the  tevel  widths  of  all  the  processes  contributing  to  the 
level  decay. 
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Assuming  the  source  lineshape,  absorber  cross-section,  and  distribution  of  nuclear 
environments  as  given  by  Eqs.  (35),  (37),  and  (36),  respectively,  the  intensity  of  the  beam 
r  absorbed  by  a  thin  absorber^  (Eo  =  0  without  loss  of  generality)  is  given  by 


j 


N. 


2,. 2  ,2, .2 

-Ej/Aj  -Ej/Aj 

e  e  Ae.  AEj 


r(r.  A,  S)  = 


I  o 

o  o 
4jcAjA2 


N, 


-.1 


X 


(r/2ji) 


r  r 
1^1 


(E-e.-Sf  +  (jrp.y  (E-ep"  +  (r/2)^ 


(38) 


where  Ns,  Na  are  the  total  number  of  sites  with  energy  shifts  Ci  and  ej  in  the  source  and 
absorber,  respectively,  and  where  na  is  the  number  of  nuclei  per  cm^  in  the  absorber. 
Converting  the  sum  to  an  integral,  one  has  (see  Ref.  24). 

r(r,A,S)  =  IoCTo(r.,/A)  n.  Re  (6^)  erfc(Z)  (39) 


where 


A^s(Aj  +  A^/2  ,and  Zs(r  +  iS)/y2  A  . 


The  maximum  normalized  resonant  absorption  (S  =  0)  is  given  by 

r(r^)  •  rcTAS = oy  (ii„n.o„) 


(40) 


r  ^ 

-^11 

~  2  {TJ 


xe*  erfex 


]. 


where 


X  = 


72  A 


Figure  42  shows  the  normalized  plot  of  expression  in  Eq.  (40). 


A  thin  MOssbauer  absorber  is  one  for  which  the  lineshape  of  the  beam  passing  through  the  absorber  is 
not  sufficiently  modified  to  show  saturation  effects.  This  is  discussed  in  ^qiendix  B,  Ref.  24,  where  it 
is  shown  that  for  the  MOssbauer  thickness  b  <  10  the  ratios  of  the  intensities  of  the  lines  in  the 
spectrum  are  close  to  their  theoretical  values  and  saturatkm  dfects  are  small.  Under  such  conditions,  it 
is  appropriate  to  use  Eq.  (34)  for  the  absorption  calculations  instead  ctf  the  general  expression  given  in 
Eq.  (41),  Ref.  24. 


NORMALIZED  RESONANCE  ABSORPTION  INTENSITY,  l"(r,A) 
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INHOMOGENEOUS  BROADENING  PARAMETER  a 


FIgura  42.  Plots  of  ths  Maximum  Normalized  Resonance  Absorption  Curves, 
r(r,A),  Calculated  From  Eq.  (40)  and  Plotted  as  s  Function  of  A/r.  The 
numbers  labeling  the  curves  give  the  different  values  of  ryr.  The 
dashed  horizontal  arrow  shows  the  effect  on  the  resonance  expected 
from  RF  pulsing  techniques  and  the  solid  arrow  the  reduction  In 
Inhomogeneous  broadening  expected  from  good  crystal 
preparation  techniques  from  ordinary  isomeric  samples 
with  a  broadening  parameter  of  a  s  10*. 
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Note  that  there  is  a  dependence  on  x  =  through  the  exponential  and  error 

function  terms  and  a  dependence  on  F  alone  through  the  factor  Fy/F.  When  the 
inhomogeneous  broadening  is  reduced  by  some  means  A  — >  0,  x  — > ««»,  and  I"  t^proaches 
a  mgyimiuTi  1/2  loOona  (Fy/F).  When  the  lines  are  homogeneously  broadened  to  obtain 
overltq)  F  ->  A  or  in  the  limit  F  x  -» «»  and  the  second  term  becomes  maximum  but 
then  (Fy/F)  0.  Thus,  although  total  overlap  is  achieved,  the  resonant  effect  disappears 
because  the  fnaximum  radiative  cross-section  on  resonance  is  reduced  to  keep  the  total 
integrated  cross-section  constant 


This  does  not  happen  when  Fy/F  =  1,  but  this  is  highly  unlikely  in  nuclear  systems. 
Generally,  the  total  width  F  is  a  sum  of  several  parts  including  the  radiative  width  Fy,  the 
internal  conversion  Fa,  and  other  transitions  Fi  so  that 


i 

Under  these  more  realistic  conditions  homogeneous  broadening  alone  will  not 
increase  the  resonant  effect  even  though  axnplete  overlap  of  lines  is  achieved  unless  we  go 
to  the  "motionally  "  narrowed  limit  and  all  the  inhomogenetnis  effects  are  wiped  out.  (See 
Ref.  24,  Fig.  24.) 


C.  COHERENT  APPROACH 

To  treat  the  problem  coherently,  since  we  do  not  know  which  nucleus  emitted  and 
which  is  going  to  be  excited,  we  have  to  sum  over  the  amplitudes  for  all  nuclear  processes. 
The  an^litudes  have  to  be  calculated  for  a  scattering  process  since  we  also  have  to  include 
the  deexcitations  coherently. 

Consider  the  emission  of  a  photon  of  energy  E  by  the  i*^  nucleus  in  the  source  at  a 
restmance  energy,  Eq  +  Ci,  where  e,  is  the  energy  shift  due  to  the  environment  This  shift 
is  responsible  for  the  inhomogeneous  broadening  of  the  system  linewidth.  The  photon 
Y  *,  scattered  by  the  nucleus  in  the  absorber,  has  energy  E'.  The  process  is  depicted  in 
Rg.  43. 

Heitler  (Ref.  25)  has  shown  how  to  calculate  the  differential  scattering  cross-section 
for  a  single  resonator  (nucleus  or  atom)  using  sectmd-order  perturbation  theory.  Boyle  and 
Hall  (Ref.  17)  have  generalized  these  results  to  a  system  of  many  nuclei. 


Our  discussion  follows  these  two  developments.  Following  Heitler  and  Boyle  and 
Hall,  the  amplitude  for  the  process  depicted  in  Fig.  43  where  the  i***  nucleus  in  the  source 
emits  a  photon  of  energy  E  and  the  nucleus  in  the  absorber  absrabs  the  photon  of  energy 


E  and  emits  a  j^otm  of  energy  E'  is  given  by 


where  Pi  is  the  probability  nucleus  i  emits  at  energy  E  and  all  the  energy  independent 
terms  have  been  lumped  into  So-  For  identical  source  and  scatterer  nuclear  environments 
(Ei  =  Ej  s  0)  Eq.  (41)  reduces  to 


■a  — — 

(E--E)-i||  (E'-Eo)-if 


2  Y  1^1  ^tI 

From  Boyle  and  Hall's  result  we  find  that  Sg  = -r - 3“  -yr  F(0)  where  y  is  the  width 

J  J 

of  the  ground  state  due  to  the  presence  of  radiation,  F  is  the  natural  width  of  the  transition 
and  the  angular  term  can  be  written  in  terms  of  CHebsh-Gordon  coefficients  and  vector 
spherical  harmonics  as  shown  in  References  26  and  27. 

For  a  ^tem  ci  N  source  nuclei  and  M  absorber  nuclei  we  get  for  the  total  scattered 
intensi^ 

of  whidi  the  incoherent  part  is  given  by 

N.M  2 

i.  j 


and  die  interference  term  Iw 


source 


absorber 


Figure  43.  A  Graphic  Representation  of  the  Scattering  Process.  Nuclei 
In  the  source  emit  photons  that  are  scattered  by  the  absorber  and 
eventually  detected.  As  on  example,  processes  f^yand 
If  are  possible  comributors  but  not  ft  ->  /  -» tf. 


Assuming  a  normal  distribution  of  nuclear  environment  as  given  by  Eq.  (36)  and 
ccmverdng  the  sums  to  integrals  we  get: 

■+mm  -H* 

Iu>ui  ®  =  J  ^  J  J  ^  E\  e,)  A* 


X  (E.  E*.  e^)  P(e)  P(e,)  PCEj)  n,n5  ^dA 


(46) 


with  die  incoming  beam  profile  given  by 
Ibeam  (E,  E)  = 


!*k 


C-'.-f  .[f ) 


A(E,  E.  e’)  =  S- 


1 


[(E-E)-i^]  [[E--e’-Eo)-ir/2] 
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(47) 


(48) 


and  ng,  na  arc  the  number  densities  of  nuclei  in  source  and  absorber,  respectively,  and  dfli 

is  the  solid  angle  intercepting  the  absorber  area  dA.  The  factor  ^  n,  picks  out  the  nuclei 
emitting  in  the  direction  of  the  absorber  and  intercepting  the  absorber  area  dA. 


For  maximum  effect  S  =  Eq  =  0.  Assuming  Eq  =  0  and  replacing  sj  by  its 
definition  we  get 


If  we  look  at  Eq,  (46a)  and  group  o©  r2/4,  y/2,1,  and  with  the  resonant 

denominators  containing  F  and  y,  respectively,  then  these  terms  wiU  be  normalized  like  an 
absorption  cross-section  with  maximum  o©  in  units  of  area  and  a  beam  with  integrated 
intensity  of  one  photon/area  [compare  Eq,  (35)  and  Eq.  (37)).  nie  maximum  intensity  will 
occur  at  energies  when  all  terms  are  on  resonance.  The  factor  multiplying  the  overall 

expression  will  be  ^  In  general  when  the  total  cross-section  F  is  increased  to  A  to  obtain 

F, 

overlap  of  resonance  lines,  the  factory  will  reduce  the  overall  effect 

We  first  perform  the  integral  over  E  in  Eq.  (46a)  using  the  result  of  Appendix  F 
(bdow)that 
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Equatirai  (46a)  can  be  written  as 


2itUjorvWAl 
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^  f  _l_rL_  _i_ 
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(46b) 


The  right-hand  side  of  Eq.  (46b)  can  be  evaluated  further  by  noting  that  the 
expression  breaks  up  into  a  product  of  three  integrals  of  form 


Bj  (z.  A,  n  =  f - ^  2"  ^ — 2  =  ^  crfc  z 

A  (E'-x)  +(rny  ^ 

B  (z  A,n=f — -iie''erfcz, 
~  (E’-x)-i[7] 

B3(z2,  A,r)  =  i crfcz2-2jtie  ^ 

iE'-T  _(iE'  +  £l 

z,r±ffl'.,. — ^  ,™iz.=J — 2.J  ^ 


erfcz 


where 


andz2' 


The  expression  for  B  was  obtained  earlier  [sec  Eq.  (39)  and  Ref.  24].  The 
expressions  for  B2  and  B3  are  derived  in  Appendix  D  (below).  Using  expressions  (49), 
(SO),  and  (51)  we  can  rewrite  (46a)  as 

I|oi(E  )  =  sj  n.nj  B,  (z.  A,  H  Bj  (z  A.  O  B,  (z,,  A,  F)  . 

»fA 

The  total  number  of  events  at  the  detector  is  obtained  by  integrating  over  the  scattered 
energy  E'  to  get 


ffl-lO 


dE’  (z,  A.  n  (2.  A.  D  B3  (Zj.  A.  O  • 


The  integration  has  to  be  done  numerically  for  artntrary  parameters  but  in  the  limit  ■ 

2A 


ii)  n,  nj 


we  get  an  explicit  result  (see  Appendix  E,  below)  I  Oo  which  can  be 

interpreted  as  follows:  For  A  ^  0  and  F  finite  (i.e.,  the  natural  linewidth)  the  maximum 
effect  is  observed.  For  A  »  finite  and  F  ->  «>  (i.e.,  large  homogeneous  broadening)  the 


effect  goes  to  zero  because  of  the  factor 
calculatitm. 


(W 


in  agreement  with  the  incoherent 


D.  CONCLUSIONS 

We  have  examined  the  effect  of  inhomogeneous  broadening  on  nuclear  resonance 
absorption  in  a  Mbssbauer  experiment  As  the  broadening  increases,  the  resonance  effect 
between  source  and  absorber  decreases.  Increasing  the  homogeneous  broadening  (i.e.,  the 
broadening  of  each  resonance  line  in  a  system  by  the  same  amount  in  the  source  and 
absorber)  increases  the  overlap  of  lines  but  decreases  die  overall  Mdssbauer  effect  This  is 
in  contradiction  to  statements  made  at  a  recent  ultrashort  wavelength  laser  meeting.  This 
happens  in  both  an  incoherent  and  a  coherent  calculation.  Thus,  the  homogeneous 
broadening  in  both  the  source  and  absorber  by  itself  will  not  in  general  improve  the 
resonance  effect 
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APPENDIX  A 

A  MONTE  CARLO  SOLUTION  OF  THE 
MAXWELL-BLOCH  EQUATIONS 


The  Maxwell-Bloch  equadcms  considered  in  this  appendix  are: 
dt 


-  =  -7N4 

31  ■‘■^1  ®3l]“ 

1=+(e;,k 

3j  +  ^l  R^jj- 

dRti  F-  +  r,  +  r_  .  + 

^  -  (N,  -  N,)  ^  ^  +  ^3,  <A-I) 

+  L,  +  r,  +  r_  _ 

(N3 - N,) ^ 

ax  2^1 

+  1 

-#=8.^3, • 

If  the  initial  state  1 0  >  is  deftned  as  the  vacuum  state  for  the  electromagnetic  field  and  an 
eigenstate  for  the  operators  Ni,  N3,  N4,  with  eigenvalues  0, 0, 1,  then  initial  and  boundary 
amditions  are 

£3^  (x,  0)  1 0  >  =  ^j(x,  0)  1 0  >  =  Nj  (X,  0)  1 0  >  =  N3  (x,  0)  1 0  >  =  0 

and  (A-2) 

E;,(0.t)  =  E^,(0.t)  =  0 

The  quantities  gi  and  ^  are  both  positive.  It  then  follows  fix>m  Eq.  (A-2)  and  the 
last  two  equations  of  (A-1)  that 

(x,  0)  1 0  >  =  R^j  (x,  0)  1 0  >  =  0  .  (A-3) 
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» 


The  terms  and  4“  appearing  in  the  third  and  fourth  equations  of  (A-1)  are 
stochastic,  independent,  Gaussian  distributed  functions  with  zero  means  and  the  second 
nooments 

<^3J  (X,  t)  (x',  f)>  =  (1/N)  <N^it)>  8(x-x*)  5(t-f)  (A-4) 

=  (1/N)  e"^*  Y  8(x-x’)  8(t-t’)  . 


The  statistical  pariimeters  given  by  (A-4)  reflect  the  fact  that  and  ^  are 
Hennitian  ccmjugate  operators. 

Acording  to  the  last  two  equatitxis. 


(A-5) 


Substituting  from  (A-S)  into  the  fourth  and  fifth  equations  of  (A-1)  leads  to  the 
differential  equations 


E:.+  |MEr-g,(N,-N,)E*- 


r, + r.  +  r_  +  i  + 

-i— J - i(Ej+iME*)  +  4 


and 


^,  + 1  he:  =  g,  (N, -  N.)  E- ^ 


r,+r,+r_  ,  + 

^  *  JL(E;  +  ipEl+^  . 


(A-6) 


In  the  spirit  of  the  Haake-Reibold  approximations  the  initial-boundary  value 
problem  for  the  Maxwell-Bloch  system  of  differential  equations  (A-1),  with  die  population 
of  level  4  totally  inverted  initially,  can  be  treated  classically  (i.e.,  as  if  the  operators  were  all 
c-numbers).  For  this  purpose  the  c-number  initial  and  boundary  values 


£3^  (X,  0)  =  (X,  0)  =  Nj  (X,  0)  =  N3  (X,  0)  =  0,  (X,  0)  =  1 

E5j(0,t)=E;,(0,t)  =  0 


(A-7) 


replace  the  operator  initial  and  boundary  values  in  (A-2).  Also,  a  corresponding  classical 
ensemble  average  replaces  any  expectation  value;  e.g.,  the  radiation  intensity 

I(t)  =  <E-(l,t)E^(l,t)> 


is  interpreted  as  a  second  moment  of  the  stochastic  electric  field  satisfying  the  substituted 
classical  irutial-boundaty  value  problem. 
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The  c-number  substitutes  for  the  operator  quantities  appearing  in  (A-1)  will  all  be 
real  functions;  hence,  a  single  function  E  will  replace  the  Hermitian  conjugate  operators 
E+  and  E"  a  single  function  R  will  replace  the  Hermitian  conjugate  operators  R+  and  R” 
and  a  single  function  ^  will  replace  the  Hermitian  conjugate  operators  and  Thus, 

E»  +  ^ME,  =  g,(N3-N,)E-i^i^^  (E,  +  inE)  (A-S) 

replaces  both  equations  in  (A-6),  and  similarly  modified,  single  c-number  equations  will 
replace  the  two  equatirais  in  (A-S)  as  well  as  the  second  and  third  equatitxis  in  (A-1). 

For  numerical  calculations  the  corresponding  differaice  equations  are: 

(X,  t)  =  (X,  g  -  hyN^  (X,  g 

N3  (x,  t)  =  N3  (x,  to)  -  2hE  (x,  g  R  (x,  g  -  hr3N3  (x,  g  +  h-j^^  (x,  to) 

Nj  (x,  t)  =  Nj  (x,  g  -  2hE  (x,  to)  R  (x,  g  -  hf^Nj  (x,  g  +  hr3N3  (x,  g 

(A-9) 

E  (x,  t)  =  E  (x,  g  +  E  (Xo,  t)  ~  E  (Xo,  g-y  Aji  [E  (Xo,  t)  -  E  (Xo,  to)l  , 

+  hAgj  [N3  (Xo,  to)  -  Nj  (Xo,  gi  E  (Xo,  g  -y hAn  (rj/2)  E  (Xo,  g 
-  h  (rj/2)  [E  (x,  g  -  E  (xo,  to)] + gjhA^  (Xo,  g 

along  with 

R(x,  t)  =  (1/Agj)  [(1  +  y  nA)  E  (x,  t)  -  E(Xo,  t)]  ,  (A-10) 

which  comes  from  (A-S).  The  quantities  h  and  A  are  the  difference  equation  time  and 
distance  increments:  (t-g,  (x-xo). 

The  stochastic  noise  function  ^(x,t),  which,  in  the  form  ^(xq,  g,  appears  only  in 
the  last  equation  of  (A-9).  It  is  calculated  by  selecting  a  random  number  frcHn  a  normal 
distributicm  with  mean  zero  and  a  variance  given  by 

o2  =  (l/N)e*‘n  .  (A-11) 

which  comes  from  (A-4),  the  delta  functions  having  been  integrated  out  in  the  last  equation 
of  (A-9). 

As  indicated,  for  every  grid  point  (Xq,  to)  an  independent  randrxn  number  is 
selected,  from  which  ^  is  obtained.  The  random  number  ^  can  be  calculated  as  follows,  in 

A-S 


tenns  of  two  random  numbers,  ui  and  U2.  chosen  independently  from  a  uniform 
distribution  over  the  interval  (0,1). 

First,  the  number  ^  given  by 

=  y-  2  In  Uj  cos  (2jc  Uj)  (A- 12) 

is  a  random  number  from  a  Gaussian  distribution  with  mean  zero  and  a  variance  of  1.  This 
can  be  seen  from  the  fact  that  the  modulus  of  a  complex  number  for  which  the  real  and 
imaginary  parts  are  independent,  normally  distributed  random  variables  with  identical 
means  and  variances  is  Rayleigh  distributed  over  the  interval  (0,<»),  and  the  phase  of  the 
con(q)lex  number  is  uniformly  distributed  over  the  interval  (0J2n).  Then  ^  is  given  by 

4  =  0^  .  (A-13) 

where  a  is  determined  by  (A-1 1). 
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LISTING  OF  CODE 
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dovible  precision  nl(201)  ,n3  (201)  ,e31p(201)  ,e31]n(201) 
double  precision  el(250) ,e2(250) ,e3(250) ,e4(250) ,e5(250) ,e6(250) 
double  precision  al(250} ,a2(250) ,a3(250) ,a4(250) ,a5(250) ,a6(250) 
double  precision  nlt(201)  ,n3t(201)  ,e31pt(201)  ,e31iDt(201) 
double  precision  ell(201)  ,r31p(201)  ,r31in(201) 
double  precision  xnsl(201) ,xns2(201) 
double  precision  dn3dt(201)  ,teiBpn3  (201) 
real  abc(250) ,abcv(2S0) 

double  precision  pi,sd,t,giaa,gi&l,gin3,ff,cf,xf,gl, alpha 

double  precision  f0,scale,bta,d,h,yl,y2,ml,m2,m3 

double  precision  e,ef ,efl,ef2,ee,f ,giiilO,r,tp,tf ,btal,bta2,tjp 

doxible  precision  cbeta 

doiible  precision  tjnl,tjn2,aa,x 

double  precision  s5v,s6v,e0 

double  precision  et 

double  precision  ffff,aaa 

Integer  lseed,n,in,nf , Issw, isssv, is4v, is5w, is6w,ksel, is, iss,np, 
'  Isv ,  it  1 ,  it  0 ,  nout ,  nx ,  nt ,  nd ,  nxv ,  ntv ,  i  i ,  mS 
conmon/a/el , e2 , e3 
comDon/b/e4 , e5 , e6 
comnon/c/al , a2 , a3 
coinmon/d/a4 ,  a5 ,  a6 

open(unit«3 ,  f  ile«'xio» , status»'nev' ,  fom»^  formatted' ) 
open (unlt«16 , f ile>' fp6 ' , status- 'new' , form-' formatted' ) 

pi  «  dble(4.*atan(l.) ) 

n  -  40 
m  -  2000 
t  -  100. do 

gma  -  10. dO 
gml  -  l.d-6 
gm3  -  l.d-4 
ff  -  l.d-3 
cf  -  O.dO 
nf  -  10 
btal  -  O.dO 
bta2  -  O.dO 
scale  -  .IdO 
Iseed  -  1029384567 
issw  -  1 
aaa  -  l.d4 
Isssw  -  0 
cbeta  -  O.dO 
ls4v  -  0 
s5w  -  l.dO 
s6w  -  l.dO 
eO  -  O.OdO 
is5w  -  0 
alpha  -  O.dO 
ls6w  -  0 
ksel  -  0 
sd  -  l.dO 
is  -  100 
iss  -  1 
xf  -  l.dO 
gl  -  l.OdO 
to  -  l.dO 
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np  “  250 

Isw  «  0 

j 

tjnl  ■  2. do 
tjn2  *•  3. do 
tjp  ■  -t  +  l.dO 

A 

ito  ■  a  +  2 
itl  -  a  +  4 
nout  -  n 

write (3,*)  'iseed,sd«' ,iseed,sd 

w 

write(3,*)  'n,a  ■  '»n,a 

write(3,*)  'is,iss«  ',is,lss 

write(3,*)  'xf,t  ■ 
write (3,*)  *gaa  -  'f9aa 
write(3,*)  'gal,ga3-  ',gal,ga3 
write(3,*)  'gl,ft  ■  ',gl,ff 

• 

write(3,*)  'nf,fO  -  ',nf,fO 
write(3,*)  'cf,np  -  ',cf,np 
write(3,*)  'scale, isw  ■  scale, isw 
write(3,*)  'ito, itl  -  '  ,itO,itl 

• 

write(3,*)  'nout  ■  ',nout 

%n:ite ( 3 ,  * )  'btal , bta2 , tjp' , btal , bta2 , tjp 

if  (isssw.ne.O)  write(3,*)  'isssw,cbeta  ■  ',isssw,cbeta 

write (3,*)  'issw  ■  ',issw 

if (issw.gt.O)  write(3,*)  'tjnl,tjn2,aa  «  ' ,tjnl,tjn2,aa 
if  (isAw.ne.O)  write(3,*)  'is4w  -  ',is4w 
if  (isSw.ne.O)  write(3,*)  'is5w,  alpha  »', is5w, alpha 
cf  -  cf/2.d0 

• 

tnrlte(3,*)  's5w,s6w,eO  «  ',s5w,s6w,e0 
write (3,*)  'is6w,ksel  •  ',is6w,]csel 
d  =•  t/dble  ( float  (a) ) 
h  «  xf/dble ( float ( n) ) 

nx  -  np 
nt  ■  1 
dt  »  .1 
to  -  0. 

• 

nd  *■  a/np 
do  11  l-l,np 

11 

abc(i)  «  sngl(t*float(i-l)/np) 
write (16, 997)  nx 
write (16, 997)  nt 

• 

997 

foraat(lx,i4) 
write (16, 998)  dt 
write (16, 998)  to 

998 

foraat(lx, flO.4) 
wrlte(16,996)  (abc(l),  i»l,np) 

996 

f oraat ( lx, f 13 . 6) 

nxv  -  n 

ntv  -  a/is 

dtv  ■  t/ float (ntv) 

tOv  ■  0. 

do  111  i«l,n 

• 

111 

abcv(i)  -  float (l-l)*sngl(h) 

A 

eO  ■■  eO*d 
do  1000  nn«l,nf 

do  1  j«l,n+l 
nl(j)  -  O.dO 
n3(j)  -  O.dO 
n3t(j)  -  O.dO 
nlt(j)  -  O.dO 
•31p(j}  «  O.dO 

• 
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• 

e31pt(j)  -  O.dO 
«31]Bt(j)  •  O.dO 
•31n(j)  •  O.dO 

If  (l86w.eq.2.and. j .eq.ksel)  then 
e31P(j)  ■  eO 
e31pt(j)  ■*  eO 
e3lmt(j)  "  eO 
e31n(j)  «  eO 
endif 
continue 

if  (isOw.eq.l)  then 

if  (nn.eq.l)  eO  >  eO/d 
e31p(l)  >  eO 
e31m(l)  •  eO 
e31pt(l)  -  eO 
e3lBt(l)  «  eO 
endif 
ii  -  0 

do  100  i>l,iD 

if  (isSw.ne.O)  gl=dexp{-i*d*alpha*gn3/2.dO) 
if (issw.eq.O)  then 

r  ■  dexp(-gina*float (i-1) *d/2. ) *dsqrt(gtta*ff) 
if  ( (itO.lt. i. and. i. It. itl) .and. isw.eq.l)  r  *  0. 
if  (isw.eq.2)  r  >  sqrt(gma*ff) 
endif 

if (issw.eq.l)  then 

if ( float ( i-1 ) *d . It . t j nl . or . float ( i-1 ) *d . gt . t j  n2 )  then 
r  -  dexp(-gi&a*float(i-l) *d/2. ) *dsqrt(gma*ff) 
else 

r  -  dsqrt(  (dexp(-gma*float(i-l)  *d)  +  aa)  *(gina*ff) ) 
endif 
endif 

if (is8W.eq.2)  then 

if ( float ( i-1 ) *d . It . t jnl . or . float ( i-1) *d . gt . t j n2 )  then 
ffff  ■  dexp(-gm3*float(i-l) *d) 

r  »  dsqrt(  (l.dO-dexp(-gma*float(i-l)*d)  )*(gn»3*ff)  )*ffff 
else 

r  ■  dsqrt(  ((l.dO-dexp(-gma*float(i-l)*d) )  +  aa)*(g]n3*ff) ) 
endif 
endif 

if (issw.eq.3)  then 

r  -  dsqrt(  (dexp(-gina*float(i-l)*d)*gina  + 
(l.dO-dexp(-gina*float(i-l)  *d) )  * 
aaa*gm3*dexp(-gm3*float(i-l) *d) ) *ff ) 

endif 
ee  >  O.dO 
ef  -  O.dO 
do  2  j*l,n 

if  (dble (float (i) ) *d. It. tjp)  then 
bta  •  btal 

else 

bta  >  bta2 
endif 

if  (isssv.ne.O)  bta  -  bta*dexp (-cbeta*d*dble( float ( i) ) ) 
if  (is4v.eq.l.and.i.ge.2)  then 

bta  -  bta*g&3/  (gm3-t-2  .dO*dabs  (tenpn3  (n) ) ) 
endif 

if  (is4w.eq.2.and.i.ge.2)  then 

bta  *  bta*gm3/  (gm3-»-2  .dO*dabs  (teapn3  ( j ) ) ) 
endif 


e  -  s5w*(d/h)*(e3lp(j)*(e31in(j+l)-e31iQ(j)) 

+  e31m(j)*(e3lp(j+l)-e31p(j)))/gl 

+  cf*(e31p(j)  +  e31in{j))**2*d/(2.d0*gl) 

ee  *  ee  -i-  e*e 

efl  -  -  gml*nl(j)*d  +  giii3*n3(j)*d 

ef2  -  -  gin3*n3(j  )*d  +  gina*dexp(-giBa*  (i-1)  *d)  *d*fO 

ef  -  ef  +  efl*efl  +  ef2*ef2 

nlt(j)  -  nl{j)  +  e  +  efl 

n3t(j)  -  n3(j)  -  e  +  ef2 

if  (nlt(j) .gt.l.dO)  nlt(j)  -  l.dO 

if  (n3t(j) .gt.l.dO)  n3t(j)  -  l.dO 

yl  -  dble(urand(iseed) ) 

y2  *  urand(iseed) 

rnl  -  f(yl,sd)*r 

rnl  "  dcos(y2*2.d0*pi)*rnl 

rn3  «  rnl 

xnsl(j)  ■  ml 

xns2(j)  “  ml 

gmlO  ■>  gm3*(l.d0  +  bta  +  alpha) 
ml  -  s6w*ml 
m3  -  s6w*m3 
tp  -  (l.dO  +  cf*h/2.d0) 
tB  »  (l.dO  -  cf*h/2.d0) 

e31pt(j+l)  ■  (e31pt(j)*tm  +  e31p(j+l)*tp  -  e31p(j)*tm) 

+  ((n3(j)-nl(j))*e31p(j))*d*h*gl 
-(gml0/2.d0)*(e31p(j+l)-e31p(j)+cf*h*e31p(j))*d  +  m3*d*h*gl 
e31pt(j+l)  -  e31pt(j+l)/tp 

if ( (ls6w.eg.2.or.is6w.eq.5) .and. j+l.eq.ksel)  then 
if  (is6v.eq.5)  then 

et  “  dsqrt(6.5d-10)*dsqrt(n3t(j)*h*gin3*5.d0)*d 
else 
et  ■  eO 
endif 

e3 Ipt ( j + 1 ) -el Ipt ( j +1 ) +dexp ( -gm3  *d*dble ( float ( i ) ) /2 . dO ) *et 
endif 

if (is6w.eq.3.or.is6w.eq.4)  then 
et  ■  n3t(j) *h*gm3 
if (is6w.eq.4)  et  «  -l.dO 

e31pt( j+l)»e31pt(j+l)+dexp(-gin3*d*dble(float(i)  )/2.d0) *et 
endif 

e31nt(j'«'l)  «  (e31iBt(j)*tm  +  e31iB( j+1) *tp  -  e31m(j)*tm) 

+  ((n3(j)-nl(j))*e3lBi(j))*d*h*gl 
-(gml0/2.d0)*(e31in(j+l)-e31ia(j)+cf*h*e31in(j)  )*d  +  ml*d*h*gl 
e31int(j+l)  ■  e31int( j+l)/tp 

if ( (is6w.eq.2.or.is6w.eq.5) .and. j+l.eq.ksel)  then 
if  (is6w.eq.5)  then 

et  *  dsqrt(6.5d-10)*dsqrt(n3t(j)*h*gin3*5.d0)  *d 
else 
et  ■■  eO 
endif 

e3  Imt  ( j +1 )  “el  Int  ( j +1 )  +dexp  ( -gin3 *d*dble  ( float  ( i ) )  /2 .  dO )  *et 
endif 

if (is6w.eq.3.or.is6w.eq.4)  then 
et  «  n3t(j)*h*gin3 
if (is6v.eq.4)  et  «  -l.dO 

e31int(j+l)»e31int(j+l)+dexp(-gin3*d*dble(float(i)  )/2.d0)  *et 
endif 
continue 

•31p(l)  -  e31pt(l) 
e3lB(l)  *•  e3lBt(l) 


3 


333 


100 


30 


31 

1000 


35 


999 

97 


if  ((1/500) *500. eg. i)  write(6,*)  nn,l,ee,ef 

do  3  j-l,n 

nl(j)  -  nlt(j) 

n3(j)  -  n3t(j) 

e31p(j+l)  ■  e3lpt(j+l) 

e3ln(j4-i)  -  e31mt(j-t-l) 

teiBpn3(j)  -  e31pt(j)  *(e3lpt(j+l)  -  e31pt(j))/h 
eii(j}  ■  e31pt( j) *e3lmt(j) 
do  333 

r31p(j)  -  (e31p(j+l)  -  e3lp(j))/h 
r31ffl(j)  -  (e31in(j+l)  -  e3lin(j))/h 
if  ( ^/is) *is.eg.i)  then 
e31p(l)  *  scale  -t-  e31p(l) 
eil(l)  «  scale  ->■  eil(l) 
r31p(l)  *  scale  +  r31p(l) 
e31p(l)  -  -scale  +  e31p(l) 
eli(l)  «  -scale  +  eii(l) 
r31p(l)  »  -scale  +  r31p(l) 
endlf 

if  ( ( (i-l)/nd) *nd.eq.i-l)  then 
el (11+1)  ■  nl(nout-l) 
e2(li+l)  «  e31p(nout) 
e3(il+l)  ■«  n3(nout-l) 
e4(il+l)  ■  fO*dexp(-gina*float(i)  *d) 
e5(li+l)  -  r31p(nout-2) 
e6(ii+l)  ■  e3lp(nout)*e31m(nout) 
ii  -  ii  +  1 
endif 

if  (Biod(i,10)  .eq.0)%n:lte(*,*)i 
continue 

if  (nn.eq.l)  then 
do  30  i«l,np 
al(i)  -  el(i) 

a2(i)  -  e2(i) 

a3(l)  -  e3(i) 

a4(i)  -  e4(i) 

a5(i)  -  e5(i) 

a6(l)  -  e6(i) 

else 

do  31  1*1, np 
al(i)  -  al(i)  +  el(i) 

a2(i)  -  a2(i)  +  e2(i) 

a3(i)  -  a3(i)  +  e3(i) 

a4(l)  «  a4(l)  +  e4(i) 

a5(i)  -  a5(i)  +  e5(i) 

a6(l)  -  a6(l)  +  e6(i) 

endif 
continue 
do  35  1*1, np 
al(i)  -  al (i) /float (nf) 
a2(i)  -  a2 (i) /float (nf) 
a3(l)  -  a3 (i) /float (nf) 
a4(i)  «  a4 (1) /float (nf) 
a5(i)  -  a5 (i) /float (nf) 
a6(i)  -  ae (i) /float (nf) 
mS  ■  (m+l)/ls 

write(l6,999)  (a6(i),  i-l,np) 
write(16,999)  (a6(l),  i-l,np) 
f  ormat  ( lx ,  f  2  0 . 15 ) 
forBat(2el5.6) 
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99 


format(5el3.5) 

X  -  O.dO 
do  777  i*l,np 
777  X  “  X  +  a6(i) 

X  *  x*t/dble ( float (np) ) 

write(3,*)  'integral  of  intensity  »  ',x 

write(6,*)  'iseed,  int  of  int  -  iseed,x 

stop 

end 

doiible  precision  function  f(y,s) 
double  precision  y,s 
f  -  s*dsqrt(2.*dlog(l.d0/y) ) 
return 
end 

FUNCTION  URAND  (lY) 

INTEGER  lY 
C 

C*************************** ******************************** 

c 

C  ROBERT  RENKA 

C  NORTH  TEXAS  STATE  UNIV. 

C  (817)  565-2767 

C 

C  THIS  FUNCTION  RETURNS  A  UNIFORMLY  DISTRIBUTED  RANDOM 
C  NUMBER  IN  THE  INTERVAL  [0,1).  THE  INTEGER  SEED  lY  IS  BOTH 
C  AN  INPUT  AND  OUTPUT  PARAMETER.  IT  SHOULD  BE  INITIALIZED 
C  TO  AN  ARBITRARY  INTEGER  PRIOR  TO  THE  FIRST  CALL,  AND 
C  SHOULD  NOT  BE  ALTERED  BETWEEN  SUBSEQUENT  CALLS  UNLESS  THE 
C  SEQUENCE  OF  RANDOM  NUMBERS  IS  TO  BE  REPEATED. 

C  URAND  IS  A  MODIFICATION  OF  THE  CODE  IN  FORSYTHE,  MAL- 
C  COUf,  AND  MOLER,  AND  IS  INTENDED  ONLY  FOR  EXECUTION  ON  THE 
C  IBM  PC  OR  A  COMPATIBLE  MACHINE.  THE  ORIGINAL  VERSION 
C  FAILS  WHEN  COMPILED  WITH  THE  MICROSOFT  VERSION  3.2  COMPI- 
C  LER  BECAUSE  OVERFLOW  IN  INTEGER  MULTIPLICATION  IS  TREATED 
C  AS  AN  EXCEPTION  AND  RESULTS  IN  THE  (HEX)  VALUE  80000000 
C  RATHER  THAN  THE  TRUE  PRODUCT  MODULO  2**32. 

C 

C  INTRINSIC  FUNCTIONS  CALLED  BY  URAND  —  DBLE,  FLOAT,  MOD 
C 

c*********************************************************** 

c 

INTEGER  lAH,  lAL,  IC,  lYH,  lYL,  lYl,  IY2,  K,  KD2,  MIC, 

.  M2 

REAL  S 
C 

DATA  IAH/12867/,  IAL/63149/,  IC/453816693/,  K/65536/, 
KD2/32768/,  MIC/1693666955/,  M2/1073741824/ 

C 

C  A  SEQUENCE  OF  RANDOM  INTEGERS  IS  GENERATED  BY  lY  « 

C  (IA*IY  -I-  IC)  MODULO  M,  WHERE  lA  AND  IC  ARE  CHOSEN  TO 
C  SATISFY  THE  PROPERTIES  RECOMMENDED  BY  KNITTH  (VOL.  2) , 

C  AND  M  -  2**31  IS  THE  LARGEST  REPRESENTABLE  INTEGER 
C  (PLUS  1)  .  lAL  AND  lAH  ARE  THE  LOW  AND  HIGH  ORDER  HALVES 
C  OF  lA,  M2  -  M/2,  AND  MIC  -  M  -  IC. 

C  COMPUTE  A  SCALE  FACTOR  S  FOR  CONVERTING  lY  TO  TYPE  REAL. 

C  S  IS  THE  CLOSEST  FLOATING  POINT  NUMBER  TO  1/M. 

C 

S  -  .5D0/DBLE(M2) 

C 

C  SET  lY  TO  IA*IY  MOD  M  (OR  IA*iy  MOD  M  -  M) ,  AVOIDING 


ooooooo  oo 


OVERFLOW.  lYl  IS  IN  AND  IY2  IS  IN  [0,M). 

lYH  ■  lY/K 

lYL  ■  lY  -  K*IYH 

IF  (lYL  .EQ.  0)  THEN 

lY  -  K*MOD(IAL*IYH,KD2) 

GO  TO  1 
ENDIF 

lYl  -  2*(IYL*31574-M2)  +  lYL 

IY2  -  KD2*MOD(MOD(IAH*IYL,K)  +  2*MOD(IAL*IYH,KD2) ,K) 
IF  (IYl/2  +  IY2/2  .LT.  M2)  THEN 
lY  -  lYl  +  IY2 

IF  (lY  .LT.  0)  lY  -  (lY  +  M2)  +  M2 
ELSE 

lY  -  (lYl  -  M2)  +  (IY2  -  M2) 

ENOIF 

THE  FOLLOWING  STATEMENT  IS  DESIGNED  TO  AVOID  OVERFLOW  IN 
THE  FOLLOWING  INTEGER  ADDITION. 

1  IF  (lY  .GE.  MIC)  lY  -  (lY  -  M2)  -  M2 

SET  lY  TO  (lY+IC)  MOD  M  AND  MAP  lY  TO  THE  UNIT  INTERVAL. 

lY  -  lY  +  IC 

URAND  -  FLOAT (IY)*S 

RETURN 


NOTATION  FOR  MONTE  CARLO  CODE  ST18 


In  Paper 

In  Code 

Interpretation 

N1 

nl 

normalized  population  of  lower  level 

N3 

n2 

normalized  population  of  upper  level 

N4 

exp  (-It) 

normalized  population  of  storage  level 

r31p 

polarization 

e; 

c31p 

electric  field 

Es'i 

r31m 

polarization 

c31in 

electric  fiekl 

Y 

gma 

punq[)ing  rate 

n 

gtnl 

decay  rate  of  lower  level 

r3 

gm3 

decay  rate  of  upper  level 

N 

m 

number  of  cooperating  nuclei 

gl 

gl 

coupling  constant 

cf 

attenuaticxi  coefficient 

ml,  xnsl 

noise  source 

^31 

m2,  xns2 

noise  source 

t 

total  time  intnval 

xf 

length  of  space  interval 

n 

number  of  space  steps 

m 

number  of  time  steps 
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In  Code 


p,  np 


number  of  points  in  graphs 


r<p=ar3 

r0=br3 


alpha 


^nl,  tjn2 
cbeta 


isssw 


standard  deviation  in  noise  generator 

seed  for  random  noise  generation 

time  steps  skipped  between  i/o 

issw  =  0,  pump  noise  only  (Haake-Reibold),  Fi(t), 

issw  =  1,  regular  pump  noise  plus  additional  non 
decaying  term,  Fi(t)  +  Gi(t) 

issw  =  2,  spontaneous  decay  noise  only,  F2(t) 

issw  =  3,  pump  noise  plus  spontaneous  decay  noise, 
Fi(t)  +  F2(t) 

inhomogeneous  broadening  parameter 
homogeneous  broadening  parameter 
b  value  on  (0,  ^p) 
b  value  on  (#,t) 

for  issw  *  0,  term  for  modifying  noise 

if  issw  ^  0,  noise  changed  on  (tjnl,  tjn2) 

exponential  degradation  coefficient  of  btal  or  bta2 

isssw  =  0,  no  exponential  degradation  of  btal  or  bta2 

isssw  =  1  for  exponential  degradation  of  btal  ot  bta2 
value  on  (^p,  t) 

dephasing  rate  from  inhomogeneous  broadening 
dephasing  rate  from  homogeneous  broadening 
initial  value  of  N4 
isw  =  0,  regular  noise 
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isw  =  1,  noise  shut  off  for  itO  through  itl  time  steps 
isw  =  2,  no  exptmential  noise  degradation 
nout  space  step  at  which  data  is  collected  (normally  nout  =  n) 

isSw  is5w  =  0,  no  time  dependence  on  gl 

is5w  =  1,  time  dependence  on  gl  ->  gl  e”®^ 3*/^ 
is6w  is6w  =  0,  noise  used  to  initiate  pulse, 

boundary  c(»idition  on  field  at  x  =  0  £0  =  0 

is6w  =  1,  boundary  condition  on  electric  field  at  x  =  0, 

eo^tO 

£0  boundary  c<»idition  on  electric  field,  E  =  £o^-^3^ 

sSw  s5w  =  1.0.  nonlinear  ten  in 
sSw  =  0.0,  nonlinear  ten  out 
s6w  s6w  =  1.0,  noise  active 

s6w  =  0.0,  noise  not  active 


Input  noise  is  generally  a  product  of  a  term  randomly  varying  in  time  and  a  constant  factor 
or  a  term  exponentially  varying  in  time.  The  random  term  is  obtained  from  a  random 
fimction  generator  and  the  terms  use  fix*  diffemit  values  of  ISSW  are 

I^(t)  =  (l-c"^)A^c’^^‘ 

g.(,).aX 
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INPUT  PARAMETERS  FOR  SAMPLE  RUN  (ST18  CODE) 


iseed .  sd=  10293: 
n .  m  = 
is .  iss= 

xf.t=  l.OC 

gma  =  10 .  C 

gml .  gm3=  1 .  OC 

gl .  ff  =  1 .  OC 

nf.fO  = 
cf .  np  =  0 .  (K 

scale .  isw  =  0 

itO .  itl  = 
nout  = 
btal .  bta2 .  tjp 

101 . 000000000000 
’"isssw.cbeta 
*is4w 

*is5w .  alpha 
*is6w .  eo 
issw  s 
tjnl.tjn2.aa  = 


1029384567  1.00000000000000 

40  2000 

100  1 
1 . 00000000000000 
10.0000000000000 
1 .  OOOOOOOOOOOOOOE-006 
1 . 00000000000000 
10  1 .00000000000000 
0 .  OOOOOOOOOOOOOOE-000 
0 .  lOOOOOOOOOOOOOOE 
2002  2004 

40 

0 .  OOOOOOOOOOOOOOE-000 


100.000000000000 

1 .  OOOOOOOOOOOOOOE-004 
1 .  OOOOOOOOOOOOOOE-003 


0 .  OOOOOOOOOOOOOOE-000 


0 .  OOOOOOOOOOOOOOE-000 


2.00000000000000 


3.00000000000000 


♦parameters  isssw.cbeta,  is4w,  is5w,  alpha,  is6w,  eo  printed  out  only  if  non  zero 
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APPENDIX  B 


MEAN  INTENSITY  IN  THE  LIMIT  OF  FAST  PUMPING 


B-l 


APPENDIX  B 

MEAN  INTENSITY  IN  THE  LIMIT  OF  FAST  PUMPING 


At  early  times  the  Maxwell-Bloch  equations  are  effectively  linear  because  of  the 
initial  ctmditions 

E+  (X,  0)  1 0  >  =  R-  (x.  0)  1 0  >  =  0,  (B-1) 

which  make  the  nonlinear  terms  quadratically  anall.  The  resulting  linearized  equations  are 
solvable  in  closed  fornL 

With  die  boundary  condition 

E+(0.t)  =  E-(0,t)  =  0  (B-2) 

inclosed  at  the  left  end  of  the  active  volume,  taking  into  account  the  statistical  properties  of 
the  stochastic  source,  the  mean  radiation  intensity  at  the  right  end  of  the  active  volume  is 
dien  given  by 

I(t)  =  <E  (l,t)E"(l,t)> 

t 

-  (yN)J t{I„(2[F{t,  -  {I,  (2[F(t.  ,  (8-3) 

o 

where  the  ]o,i(t)  are  modified  Bessel  functions  and 

F(t,t’)  =  t-t'  +  (c'^*-c'^‘)/y  .  (B-4) 

In  the  limit  of  fast  pumping,  y  ->  oo,  the  average  intensity  given  by  (B-3)  becomes 

L  (t)  =  (l/N){[Ioa/t)]^-[Ij(2/I)]^}  .  (8-5) 

For  the  time  rate  of  change  W(t)  of  the  probability  of  finding  the  mean  intensity  I(t) 
above  a  reference  level  I^ef,  Haake  and  Reibold  give  the  equation 

W(t)  =  ti(t)I^(t)^  exp  [-1^(0]  .  (B^) 
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Haake  and  Reibold  remark  that  W(t)  should  be  a  reasonable  approximation  to  the 
probability  density  of  the  first  passage  time  of  I(t)  above  and  therefore  at  least  a  rough 
approximation  to  the  delay  time  probability  density. 

The  time  at  which  the  maximum  value  of  W(t)  occurs  is  then  an  estimate  of  the  most 
probable  delay  time  of  the  radiated  pulse.  Some  typical  delay  time  estimates  based  on  B-6) 
with  a  value  of  1  assigned  to  and  various  values  of  N  are  given  in  Figs.  (B-l)-(B-6), 
which  illustrate  the  conesponding  function  W(t). 
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W(t) 


Delay  s  9.4 
Peak  =0.20 


Figure  B*1.  Passage  Time  Probability  Density  In  the 
Limit  of  Infinitely  Fast  Pumping:  N  s  1000 


t  :  =  15.15.5  ..30 

0.2 


W(t) 


Delay  =  194 
Peak  B  0.15 


Figure  B*2.  Passage  Time  Probability  Density  in  the 
Limit  of  infiniteiy  Fast  Pumping:  N  s  10^ 


- ! 


HM-i 


W(t) 


Till 


I  I  I 


J  Delay  =  32.3 

J  Peak  s  0.12 


Figure  B*3.  Passage  Time  Probability  Density  in  the 
Limit  of  Infinitely  Fast  Pumping:  N  s  10^ 
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Figure  B<4.  Passage  Time  Probability  Density  In  the 
Limit  of  infinitely  Fast  Pumping:  N  &  tO^ 


Figure  B-5.  Passage  Time  Probability  Density  In  the 
Limit  of  Infinitely  Fast  Pumping:  N  =  lO'''' 


t 


90 


Delay  e  88.3 
Peak  s  0.07 


Figure  B-6.  Passage  Time  Probability  Density  In  the 
Limit  of  Infinitely  Fast  Pumping:  N  = 
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APPENDIX  C 

AVERAGE  OF  A  FUNCTION  WITH  A 
STOCHASTIC  ARGUMENT 


This  appendix  details  an  investigation  of  the  following  problem:  to  calculate  the 
stochastic  average  of  the  functional 


where  g(t)  and  h(t)  are  identical  stochastic  functions,  each  of  which,  as  a  stationary  Maikov 
process,  independently  takes  on  random  values  ±  in  successive  infinitesimal  tinyi  intervals 
At  The  calculation  is  equivalent  to  finding  the  avmge 


where  f(t)  is  a  stochastic  functitm  of  the  sanoe  type  as  g(t)  and  h(t),  but  which  can  take  on 
any  of  die  values 

ai  =a-p,a2  =  a  +  p,a3  =  -a2,a4=-ai  ,  (C-1) 


with  equal  probability. 

As  in  the  case  treated  by  Ref.  C-1, 


lj>  . 


(C-2) 


where  F  is  a  4*b  order  diagonal  matrix  with  the  main  diagonal  elements  ai,  given  by  (1) 
and  W  is  the  4^  order  matrix 


C-3 


(C-3) 


-2a 

Cl, 

Cl, 

o" 

w  = 

a. 

-2Cl, 

0, 

Cl 

a. 

0, 

-2Cl, 

Cl 

-  0, 

Cl, 

a-2fij 

According  to  Ref.  C-2,  each  off-diagonal  element  Wij  of  W  is  the  conditional  probability 
per  unit  time  of  f(t)  taking  the  value  aj,  at  time  t  +  At,  given  that  it  had  the  value  ai  at 
timet 


The  zero  elements  occur  because  the  corresponding  conditional  probabilities  are 
second  order  in  At  The  diagonal  elements  all  have  the  value  -2Cl  because,  as  Ref.  C-2 
observes,  the  fact  that  for  each  condition  the  sum  of  conditional  probabilities  must  be  one 
inq>lies  that  the  sum  of  the  corresponding  probability  rates  must  be  zero. 

The  straightforward  way  to  evaluate  (C-2)  is  to  diagonalize  the  matrix  iF  +  W.  The 
first  step  is  to  calculate  its  eigenvalues  Xi  by  solving  the  equatim 

det(iF  +  W-I)  =  0  ,  (C-4) 

where  I  is  the  4*  order  identity  mattix. 


As  is  evident  from  (C-3)  and  the  definition  of  F,  (C-4)  is  equivalent  to 


det 


iaj-4,  Q,  Q,  0  ” 

iOj-H,  0,  Q 

Q.  0,  iOj-ji,  Q 
0.  Q,  ia^-p. 


(C-5) 


where 

P  =  X  +  2  .  (C-6) 

Using  (C-1),  evaluating  the  determinant  in  (C-S)  leads  to  the  quartic  equation  in 

M.'‘  +  2(oUp^-2n^n2  +  (a»-pV  =  0  ,  (C-7) 


which  is  also  a  quadratic  equation  in  and  therefore  explicitly  solvable.  All  solutions  of 
(C-7)  are  given  by 

\i  =  ±[yfc^  ±  ;  (C-g) 

C-4 


therefore  according  to  (C-6)  the  required  eigenvalues  are  the  four  quantities 

X  =  -2Qj  [Va^-a^  .  (C-9) 

The  next  step  is  to  calculate  the  eigenvector  corresponding  to  each  of  the 
eigenvalues  given  by  (C-9).  The  components  Tik  of  the  eigenvector  corresponding  to  the 
k*h  eigenvalue  will  satisfy  the  system  of  equations 


(ioij  —  2f2)  Tjjf  +  Q  Tqv 

+  QT3k 

f2T}][  +  (iot2  “  2f2)  T2k 

+  QT4k 

=  ^T2k 

flTik 

+  (i(X3  —  2f2)  T3k 

+QT4k 

=  ^T3k 

fl[T2k 

+  QT31C 

+(io^  —  2  0)T4|^ 

=  X|f  T4j[ . 

(C-10) 

Subtracting  the  first  equatim  in  (C-10)  fix>m  the  last  and  using  (C-6)  leads  to 

=  .  (C-11) 

From  the  third  equation  it  follows  that 

T3k  =  -^(T4k  +  Tn)/(ia3-p)  ,  (C.12) 

and  from  the  first  that 

T2k  =  -  T3k  -  (»«!  -  V^)  Tik  /«  .  (C-13) 

For  each  value  of  k  the  Eqs.  (C-1 1)-(C-13)  give  T2k,  Tj^,  and  T4k  in  terms  of  Ti^.  Those 
assignments  and  the  assignments  of  any  arbitrary  complex  numbers  for  the  T^  will 
provide  solutions  for  the  system  of  Eqs.  (C-10). 

The  Tn  satisfying  (C-10)  are  elements  of  a  diagonalizing  transformation  matrix  T, 
such  that 

exp  ((iF  +  W)  t]  =  Texp(At)  T"*  ,  (C- 14) 

where  A  is  the  diagcxial  matrix  for  which  the  main  diagonal  elements  are  the  eigmvalues  Xj. 

As  just  noted,  for  each  k  the  solution  vector  of  (C-10)  will  depend  on  an  arbitrary  choice  of 
T|^.  If  this  number  is  chosen  so  that 

4 

•  <C-15) 

i=l 


C-5 


the  inverse  matrix  T"1  of  the  corresponding  T  wiU  be  equal  to  the  adjoint  matrix  "H. 


It  then  follows  that  (C-2)  is  equivalent  to 

4 

,V  »•  i>  k=l 

Since  the  a  priori  probabilities  px  of  the  four  possible  f(t)  values  are  all  equal,  each 
Pi  =  1/4  in  (C-16). 

The  relations  (C-9),  (C-11)-(C-13),  and  (C-15)  provide  the  other  quantities  in 
(C-16).  In  the  special  case  £2  =  0,  these  lelaticms  lead  immediately  to  the  result 

G  (a,  p,  t)  =  cos  (at)  cos  (pt)  ,  (C-17) 


t 

exp 

ijf(t')  df 

-  o 

where 


[ 

t  1 

G(a,  P,t)  =  ^ 

l^exp 

ij  f(t')dt'  ► 

_  o 

is  given  by  (C-16). 

The  functitm  R(t)  is  an  average  reduction  factor  defined  by 

R(t)  =  Kj  jG(a,p,t)exp[-(a2+p^]/o2]ddP  , 

O  O 

where  K  is  chosen  so  that 

R(0)=1  . 

For  the  special  case  £2  =  0,  it  follows  from  (C-17)  that 

R(t)  =  exp  (- a2t2/4)  . 

Figures  C-1  to  C-10  show  R(t)  for  various  values  of  £2  and  a. 


(C-18) 


(C-19) 


C-6 


1 


Figure  C-3.  Average  Reduction  Factor  as  a  Function  of  Time  After  Inversion. 

a  s  0.1,  o  s  0.1. 


FIgura  C^.  Average  Reduction  Factor  as  a  Function  of  Time  After  Inversion. 

a  s  100,  <T  s  10. 
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APPENDIX  D 

NOTE  ON  A  GAUSSIAN-LORENTZIAN  INTEGRAL 


The  W(z)  function  defined  in  Ref.  D- 1  can  be  used  to  evaluate  the  integral 

f  e 

I(E,r)  =  - £ -  dx  .  (D-1) 

1  (x-E)-iL 

The  W(z)  function  is  defined  by 

W(z)  =  e~^  erfc  (-iz),  Im(z)  >  0  05-2) 


and  has  the  integral  represottation 


W(z)  = 


dt 


0>3) 


for  z  in  the  upper  half  plane,  so  that 

I(E  JO  =  i  Jt  e~^  erfc  (-iz)  , 


for 


z  = 


(I>4) 


In  (D-4)  the  definition  of  z  insures  that  it  is  in  the  upper  half  of  the  complex  plane, 
as  (D-3)  requires.  However,  for  the  integral  associated  with  I(E,-r),  z  is  in  the  lower 
haliplane,  which  violates  that  condition. 

On  the  other  hand,  deforming  the  integration  contour  in  (D-3)  leads  to  the 
conclusion  that 

I(E,-r)  =  i  jc  ^  (z)  ,  (D-5) 


D-3 


where 


(D-6) 


as  long  as 


e<Im(z)  . 


But  finom  (D-3)  and  (D-6)  it  is  clear  diat 

i  jc  W(z)  -  i  Jt  ^(z)  =  ^  dx  =  2  tc  i  e“^  . 


Combined  with  (D-2)  and  (D-S).  (D-7)  inches  that 

,2 

I(E,  -  n  =  i  c  erfc  (-iz)  -  2  it  i  e 


(D-7) 
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APPENDIX  E 

INTEGRAL  EVALUATION  I 


In  the  limit 


the  expressitm 


r  oo  or  A  ->  0 


I  =  icj^2Aj  -^J  crfc  2  •  e  '  crfc  Zj  e**  erfc  Zj  -  2jcie  ^  dE  .  (E-1) 


where 


iE-f 


lE  +  y 


r  •  iE  2 

*'  = — r-^.z,  = - r 

2A  1  A  ■*  A 


(E-2) 


can  be  evaluated  explicitly. 

The  first  step  is  to  make  the  change  of  variable  in  the  integral  (xi  the  right-hand  side 
of  (E-1),  after  which  the  equation  for  I  becomes 

I  =  k(^]  rjtje**  erfcz'  e**  erfc^j^e^  erfc2^-2  jcie  ^  dE' ,  (E-3) 


where 


z'  =  i(l  +  iE'),2'j=^(2iE'-l).z;j  =  -i(2iE'-i-l)  andE’  =  ^  . 


The  next  step  is  to  introduce  the  asymptotic  limit 


lim  e*^  erfc  z  =  — ^  , 
zfn 


(E-4) 


which  (cf.  Ref.  E-1)  is  valid  for 


I  argzl< 


Then  the  equatioi  for  I  becomes 

f _ dT _ 

/W  J  (l+iE')(2iE'-l)(2iE'+l)  ’ 


(E-5) 


Evaluating  the  integral  on  the  right-hand  side  of  (E-S)  by  residues  is  straightforward  and 
gives  the  result 

T_  kJk 
^  • 
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APPENDIX  F 

INTEGRAL  EVALUATION  II 


Writing 


(F-1) 


(F-2) 


after  factoring  the  denominatOTS  in  the  integrand  on  the  right-hand  side,  Eq.  (F- 1)  becomes 


_ dz _ 

[z-(aj+ibj)]  [z-(aj-ibj)]  [z-(a2+ib2)]  [z-fa^-ib^)] 


(F-3) 


Evaluating  the  integral  on  the  right-hand  side  of  Eq.  (F-3)  by  residues  after  closing 
die  real  axis  contour  in  the  upper  half  of  the  oxnplex  plane  leads  to 

/ _ 1 _ + _ I _ 

I  =  2ra  1 2ibj[(aj-a2)  +  ii\-b^)]  [(a^-aj)  +  Kb^-i-bj)]  2ib2[(a2-a,)  +  i(b2+bj)] 

_ I _ + _ I _ 

b,[(a,-a2)+i(bj-b2)]  [(aj~a2)+i(bj+b2)]  b2[(aj-a2)+i(b,-b2)]  [(aj-a2)+i(bj+b2)]  ’ 

bitCara^)  i(bl^^b2)]  +  b2[(a^-a2)  -  i(b^+b2)]  ] 

bibjlCai-az)  +  i(bj-b2)]  [(aj-a2)^  +  (bj+b2?]  | 

I  (b,-fb2)  (8^-82)  i(b^4-b2)  (b^-b2) 

”1  b,b2[(a,-a2)  +  i(bj-b2)]  [(aj-a2)^  +  (bj+b2)^  ’ 

which  reduces  to 


^c(b^^•b2) 

bjb2[aj-a2f+(bj+b2)^ 


(F-4) 


F-3 


Substituting  from  Eq.  (F-1)  into  Eq.  (F-4)  gives  the  value  of  the  integral  on  the 
right-hand  side  of  Eq.  (F-1),  i.e., 


2jt(r-t7) 


F-4 


